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Abstract 

We present a flavour-covariant formalism for scalar triplet leptogenesis, which takes 
into account the effects of the different lepton flavours in a consistent way. Our main 
finding is that flavour effects can never be neglected in scalar triplet leptogenesis, even in 
the temperature regime where all charged lepton Yukawa interactions are out of equilib¬ 
rium. This is at variance with the standard leptogenesis scenario with heavy Majorana 
neutrinos. In particular, the so-called single flavour approximation leads to predictions 
for the baryon asymmetry of the universe that can differ by a large amount from the 
flavour-covariant computation in all temperature regimes. We investigate numerically the 
impact of flavour effects and spectator processes on the generated baryon asymmetry, 
and find that the region of triplet parameter space allowed by successsful leptogenesis is 
significantly enlarged. 
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1 Introduction 


Leptogenesis [1], i.e. the generation of a lepton asymmetry through the out-of-equilibrium 
decays of heavy particles before sphaleron freeze-out, is one of the most appealing expla¬ 
nations for the origin of the observed baryon asymmetry of the universe. As such, it 
has been the subject of a large number of studies, especially in its minimal version in¬ 
volving heavy right-handed neutrinos (for a comprehensive review, see Ref. [2]). Over 
the years, many refinements have been added to the computation of the generated baryon 
asymmetry: spectator processes M, finite temperature corrections [5] , lepton flavour ef¬ 
fects PE El EE Eo] and finally attempts to provide a full quantum mechanical formulation 
of thermal leptogenesis [fflE3E3E2E3E3EIlE8]. 

By constrast, much less work has been devoted to the leptogenesis scenarios involving 
fermionic EE [201 EH |22l [23] or scalar [2ll25l[26ll271|28l[29lEDlEI] electroweak triplets 
(for a recent review of these scenarios and additional references, see Refs. [32]). The 
CP asymmetry in scalar triplet/antitriplet decays was computed for various models in 
Refs |23, [25] . A detailed quantitative study of scalar triplet leptogenesis was performed 
in the single flavour approximation in Ref. [26], and was extended to the supersymmetric 
case in Ref. m ■ Flavour effects were addressed in flavour non-covariant approaches in 
Refs. [30[ [31j, and spectator processes were also included in Ref. [3T] . In this paper, 
we extend and improve previous works on flavour-dependent scalar triplet leptogenesis 
by providing a complete set of flavour-covariant Boltzmann equations using the density 
matrix formalism [33] [6]. We find that flavour covariance is a crucial ingredient of the 
computation of the generated baryon asymmetry, and that flavour effects also matter in 
the temperature regime where charged lepton Yukawa interactions are out of equilibrium. 
We show in particular that, as opposed to the standard leptogenesis scenario with right- 
handed neutrinos, the single flavour calculation does not provide a good approximation 
to the full flavoured computation, even when charged lepton Yukawa couplings can be 
neglected. 

The paper is organized as follows. The framework and the basic ingredients of flavour- 
dependent scalar triplet leptogenesis are presented in Section [2] In Section [3] we justify 
the use of a flavour-covariant formalism and derive the Boltzmann equation for the density 
matrix in the closed time path formalism. Spectator processes and chemical equilibrium 
are discussed in Section E Section [5] provides the set of Boltzmann equations relevant to 
the different temperature regimes, both in the flavour-covariant approach with spectator 
processes included, and in various approximations neglecting flavour covariance and/or 
spectator processes. In Section [6] we study numerically the impact of flavour effects and 
spectator processes on scalar triplet leptogenesis, and we compare the flavour-covariant 
with the flavour non-covariant computations. We find that the generated baryon asym¬ 
metry can be significantly enhanced (up to several orders of magnitude) by the proper 
inclusion of flavour effects and spectator processes, thus enlarging the region of triplet 
parameter space allowed by successsful leptogenesis. Finally, we present our conclusions 
in Section E The formulae for the space-time densities of reactions used in the paper can 
be found in Appendix lAl 
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Basic ingredients of scalar triplet leptogenesis 


2.1 The framework 


We work in the framework of the type II seesaw model {34] . i.e. the Standard Model 
augmented with a massive scalar electroweak triplet which couples to left-handed leptons 
and to the Higgs boson as follows: 


C = ~ (f a pllCia 2 Alp + fiH T ia 2 A^H + h.c.) - M\ tr(A*A), 
where C is the charge conjugation matrix defined by = — 7 ^, and 

A f = / 


( A+/V2 A++ \ 

V A 0 -A+/V2 ) ’ 


A-/y/2 A 0 * 

A— -A~/y/2 


( 2 . 1 ) 


( 2 . 2 ) 


When decomposed on the components of each electroweak multiplet, this Lagrangian 
becomes 


C = (vaCv/3 A 0 - V2v^Cep A + - e^Cep A ++ ) 

+ n A 0 * -V2H°H + A~ + H + H + A"') + h.c.} 

- M\ (A°*A° + A _ A + + A A ++ ) . (2.3) 


The triplet gives a contribution to the neutrino mass matrix: 

1 v 2 

( m A.)a/3 — 2 m 2 1 


(2.4) 


where v = (H°) = 174 GeV is the Higgs boson vacuum expectation value. In addition to 
generating small neutrino masses, the heavy scalar triplet can also create a lepton asym¬ 
metry through its decays, like the right-handed neutrinos of the type I seesaw mechanism. 
There are however important differences between the standard leptogenesis scenario in¬ 
volving heavy Majorana neutrinos and scalar triplet leptogenesis. First A has 2 types of 
decays, into Higgs bosons and into antileptons, with tree-level decay rates 

r(A - «> = ik <r(fP)M A , UA -> HH ) = A. . ( 2 .5) 

It will prove convenient in the following to introduce the quantities: 

^ = a/ tr(//t), \ H = -M-, ( 2 . 6 ) 

which control the tree-level decay width and branching ratios of the scalar triplet: 

Ta = 22 ^ (A^ + A^) Ma , (2.7) 

B e = BR(A ->• II) = X l , B h = BR(A ->■ HH) = X " . (2.8) 

A t + a h A t + A tt 
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Second, in contrast to the heavy Majorana neutrinos of standard leptogenesis, the scalar 
triplet is not a self-conjugate state. Hence, an asymmetry between the triplet and an¬ 
titriplet abundances will arise and affect the dynamics of leptogenesis. Third, being 
charged under SU( 2)l x J7(l)y, the triplets and antitriplets can annihilate before decay¬ 
ing. In order to generate a large enough lepton asymmetry, decays have to happen with a 
rate higher or similar to gauge annihilations, which are typically in thermal equilibrium. 
This requirement seems to conflict with Sakharov’s third condition [35], but the fact that 
the triplet has several decay channels makes it possible for some of them to occur out of 
equilibrium and resolves the contradiction [26] . 

Finally, another important feature of scalar triplet leptogenesis, whose significance 
has been missed so far, lies in the fact that the scalar triplet couples to pairs of leptons 
from different generations rather than to a coherent superposition of lepton flavourtjj (as 
opposed to a right-handed neutrino). As a consequence, its dynamics cannot be described 
by Boltzmann equations involving a single lepton asymmetry. This is at variance with 
the standard leptogenesis scenario, in which neglecting the effects of the different lepton 
flavours is a good approximation at high temperature. For scalar triplet leptogenesis, a 
flavour-covariant formalism must be employed, as will be discussed in Section [3] 

In order for leptogenesis to account for the observed baryon asymmetry of the universe, 
a large enough asymmetry between the leptonic decay rates of triplets and antitriplets is 
needed. It is a well-known fact [36s, ;37j that this condition is not satisfied with the minimal 
particle content of the type II seesaw mechanism, since there is no CP asymmetry in A/A 
decays at the one-loop levell. For this to happen, another heavy state (or possibly several 
heavy states) with couplings to the lepton and Higgs doublets must be added to the 
model. If this additional particle is significantly heavier than A, it will not be present 
in the thermal bath at the time of leptogenesis and its effect can be parametrized by the 
effective dimension-5 operator [26] 

A £ eS = \^Y (tWH) c (H T i<j 2 3 lp) + h.c., (2.9) 

suppressed by A Ma , where A = /// = Ma' /\' H if the heavier particle is a scalar 

triplet A', and A = M if it is a right-handed neutrino or a fermionic triplet with mass M. 
The operator (I2.9jl also gives a contribution to the neutrino mass matrix: 


1 v 2 


( 2 . 10 ) 


In full generality one should also add an effective dimension-6 operator (which arises at 
tree level if the heavier particle is a scalar triplet m, but only at the one-loop level if it 
is a right-handed neutrino or a fermionic triplet): 


A'-Ccff = , (2.11) 

2 Indeed, except for very specific flavour structures of the parameters f a p, it is not possible to find a superpo¬ 
sition of lepton doublet flavours £a = c a £ a such that the scalar triplet couplings to leptons f a pA£ a lp can 
be rewritten as /A£a^a- By contrast, the couplings of a right-handed neutrino Y\ a N\i a H can be rewritten 
yN\t Nl H , with y = \ Y ia\ 2 and £ Nl = J2 a Y la £ a /y. 

3 We found a non-vanishing CP asymmetry only at the 3-loop level, and only in the flavoured regime. 
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Figure 1: Diagrams responsible for the flavour-dependent CP asymmetries e a p. 


where rj a p^ is symmetric under the exchanges a -H- (3 and 7 -H- 6. This operator induces 
a contribution to the flavour-dependent CP asymmetries in A/A decays that vanishes 
when summed over lepton flavours, hence it only affects leptogenesis when the dynamics 
of the different flavours is taken into account. However, it generally plays a subdomi¬ 
nant role because it is suppressed by an additional power of A (and possibly also by a 
loop factor) with respect to the operator (12.91) . A notable exception arises when (12.91) 
and (12.111) are generated by a heavier scalar triplet A' with couplings to lepton and Higgs 
doublets such that XhX' h Xi^Ma/Ma'- In this case, the operator (12.111) gives the 
dominant contribution to the flavour-dependent CP asymmetries in A decays, while the 
total CP asymmetry (to which only (12.91) contributes) is small. This scenario, dubbed 
purely flavoured leptogenesis (PFL) in the literature, has been studied in Refs. [HUES]. 
In this paper we shall stick to the less specific case of dominance of the operator (12.91) and 
omit the operator (12.111) . 


2.2 CP Asymmetries in A/A decays 

The CP asymmetries in the decays of the triplets and antitriplets are defined by 

T(A -> HH) - T(A -»• HH) 

Iff = 2 -—-—-, 


( 2 . 12 ) 


_ T(A -> 44 O - T(A -> ijp) 

e ap = -y——y- (l + A*/ 3 ), 

1 A + i A 


(2.13) 


where we included a factor 2 in the definitions of en and e aa for later convenience. The 
flavour-dependent CP asymmetries e a/ g come from the interference between the diagrams 
shown in Fig. [I] With the definition (12.131) . they can be expressed in terms of tua and 
mf{ as 


€a/3 


1 Ma / Im [(m* A ) a p(m H ) a p\ 

47t v z m/\ 


(2.14) 


4 We disagree with the claim m of a strong enhancement of the generated lepton asymmetry in PFL for low 
triplet masses, which can be traced back to an erroneous term in the Boltzmann equations of Ref. [30]. Indeed, 
the second term in Eq. (23) of Ref. |30] generates a lepton asymmetry in thermal equilibrium, thus violating 
the third Sakharov condition [35] . 
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where we have introduced 


rriA = \] tr(m^mA) 
The total CP asymmetry is given by 


e A — ^ ^a/3 

a,(3 




Im 


tr (nv A m-u) 


47 t v 2 ' rriA 

where the last equality follows from CPT invariance. 


= e# 


(2.15) 


(2.16) 


2.3 Washout processes 

The lepton asymmetry generated in triplet and antitriplet decays is partially washed out 
by lepton number violating processes. These include the inverse decays £ a £p —>• A and 
£ a £p —>• A, and the A L = 2 scatterings £ a £p e>- HH and £ a H £pH, mediated by s- and 

1- channel triplet exchange, respectively, as well as by the higher order operator (12.91) . In 

addition, some processes that preserve total lepton number redistribute the asymmetries 
between the different lepton flavours, namely the 2 lepton-2 lepton scatterings £ a £p -H- £^£$ 
and £ ol £ 1 (mediated by s- and t-channel triplet exchange, respectively), and the 

2- 2 scatterings involving leptons and triplets £ a A £pA, l a A -H- £p A and £ a tp -B- AA 
(mediated by s-, t- and u-channel lepton exchange). Since the asymmetries stored in 
different lepton flavours are subject to different washout rates, these flavour violating 
processes have an impact on the erasure of the total lepton number, and we will refer to 
them as washout processes, too. 

We define for later reference the space-time density of triplet and antitriplet decays: 

7 d = Y b( A *+ ijp) + 7(A o- £ a £p)] + 7(A o HH) + 7 (A o HH) , (2.17) 

a,/3 

as well as the following combinations of flavour-dependent scattering densities: 

7 ih = Y ^ HH) + 7 (£ a H £pH)\ , 

a,/3 

74 £ = Yj [ 2 7 ( 4//3 ^ £^£s) + 7 ( 4^7 £p£-s)\ , 

a,/3,7 ,5 

= Y^ [ 7 (^qA b> £pA) + 7 (£ q A -h- £p A) + 7 {£ a £/3 AA)] 

a,/3 

3 Flavour-covariant formalism 

3.1 The need for a flavour-covariant formalism 

Leptogenesis computations are often performed in the so-called single flavour approxima¬ 
tion, in which a single direction in flavour space is considered. This is a rather accurate 


(2.18) 

(2.19) 

( 2 . 20 ) 
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approximation in scenarios in which a heavy Major ana neutrino N\ decaying at high tem¬ 
perature is responsible for the whole lepton asymmetry. Indeed, the couplings of Ni can 
be rewritten as (with the contraction of SU(2)l indices omitted) 

- ^ Yi a Ni£ a H + h.c. = -ym Nl H + h.c., (3.1) 

a 


where y = \fYla |Tia | 2 and £nx = Yha ^la^a/y- When processes mediated by the charged 
lepton Yukawa couplings are out of equilibrium, i.e. when N\ decays in the high tem¬ 
perature regime (T > 10 12 GeV), the coherence of £n x is effectively preserved by all 
interactions in the thermal plasma and leptogenesis can safely be described in terms of a 
single lepton flavour. When the lepton asymmetry is generated at lower temperature, on 
the other hand, charged lepton Yukawa interactions enter equilibrium and destroy the co¬ 
herence of £]y 1 . The single flavour approximation is no longer appropriate, and the proper 
treatment involves a 3 x 3 matrix in lepton flavour space [ 6 ], the density matrix {Af) a p, 
whose diagonal entries are the asymmetries stored in each lepton doublet £ a , while the 
off-diagonal entries encode the quantum correlations between the different flavours. This 
matrix transforms as Ai —>• U*Ag U T under flavour rotations £ —>• U£, and its evolution is 
governed by a flavour-covariant Boltzmann equatioi0. 

However, in the case of leptogenesis with right-handed neutrinos discussed above, the 
density matrix formalism is only really needed at the transition between two temperature 
regimes [22], or in the case where several heavy Majorana neutrinos contribute to the 
generation and washout of the lepton asymmetry as in Refs. [ 6 } 03] HU [38] , or in resonant 
leptogenesis [25] [39]. Otherwise there is always a natural choice of basis in which the 
Boltzmann equation for {Af) a p can be substituted for a set of Boltzmann equations for 
1, 2 or 3 flavour asymmetries. Above T = 10 12 GeV, where all charged lepton Yukawa 
couplings are out of equilibrium, the appropriate flavour basis is {£nxi£± 1 i £ 1 . 2)1 where £±1 
and £±2 are two directions in flavour space perpendicular to £n 1 . In this basis, leptogenesis 
is well described in terms of the sole lepton asymmetry A( N . Below T = 10 12 GeV, tau 
Yukawa-induced processes like q^£ T —>• Irtr are in equilibrium and destroy the coherence 
between £ T and the other two lepton flavours. However, as long as the muon Yukawa 
coupling is out of equilibrium (which is the case if T > 10 5 6 * * 9 GeV), the coherence of £q = 
(Y\ e £ e + Y\^£^)/xJ\Y\ e \ 2 + |Yi m | 2 is preserved, and the dynamics of leptogenesis can be 
described in a 2-flavour approximation, in terms of the asymmetries A ^ 0 and Aj> T . Finally, 
below T = 10 9 GeV, the muon Yukawa interactions are in equilibrium as well, and flavour 
coherence is completely broken. Leptogenesis is then governed by Boltzmann equations 
for the three asymmetries Ae e , Af and A^ t . 

The case of scalar triplet leptogenesis is significantly different, because the scalar triplet 
does not couple to a single combination of lepton flavours in general. This makes the use 


5 Except for the A L = 2 scatterings l a H -o- £pH and t a £g HH mediated by the heavier Majorana 
neutrinos N 2 and A 3 , which are neglected in this discussion. 

6 This formalism has been extended to the case where several heavy Majorana neutrinos with hierarchical 

masses play a role in leptogenesis in Ref. [35], and a fully flavour-covariant formalism in which the quantum 

correlations between the heavy and light neutrino flavours are taken into account has been developed in Ref. [55] 

(see also Refs. [401 SI] for an earlier use of a density matrix for the heavy neutrino flavours, in the scenario 
where the baryon asymmetry is generated through CP-violating oscillations of the “heavy” neutrinos below the 
electroweak scale). 
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Figure 2: The B — L asymmetry generated in scalar triplet leptogenesis as a function of z — 
tua/T, computed following three different prescriptions: density matrix computation using the 
Boltzmann equations fl5.1j) - (l5.3|h yielding A b-l = tr(A a( g) (red solid curve)] naive computation 
using the Boltzmann equations (I5.8p ~ (l5.10p written in the charged lepton mass eigenstate basis, 
yielding (Ab^) 1 = A e + A M + A r (black dashed curve)] naive computation performed in the 
neutrino mass eigenstate basis, yielding (A b-l) 2 = Ai + A 2 + A 3 (black dotted curve). The 
parameters are J14 = 10 13 GeV, A h = 0.01 and, in the notation of Eqs. (16.31) and (16.4p . 
012 = 013 = 7 t/4, 023 = 0, = 7t/ 4 and 0 2 = tt/ 6, all other phases being zero. 

of a flavour-covariant formalism unavoidable as long as the quantum correlations between 
lepton flavours are not destroyed by charged lepton Yukawa interactions. One can still 
define formally a single flavour approximation by making the substitutions (A () a a —» A^, 
f a /3 —>• A/>, u a /3 —>• A re = y / tr(K/ct) and £ a p —> ca (where £ a p are the flavour-covariant 
CP asymmetries, to be defined later), but the resulting Boltzmann equations cannot be 
derived from the flavour-covariant ones by taking the limit of vanishing Yukawa couplings. 
Hence, even if the scalar triplets decay above T = 10 12 GeV, there is no guarantee that the 
single flavour calculation will give a good approximation to the flavour-covariant result. 

A less naive approach is to neglect flavour covariance but keep track of the differ¬ 
ent lepton flavours, and write classical Boltzmann equations for the three asymmetries 
A^ a . When the lepton asymmetry is generated below T = 10 9 GeV, this is fully justified 
because quantum correlations between the different lepton flavours are destroyed by fast 
interactions induced by the muon and tau Yukawa couplings. In this case, the appropriate 
flavour basis is the charged lepton mass eigenstate basis: A e a = {Ae e ,Ag,Ag T }. When 
all charged lepton Yukawa interactions are out of equilibrium, however, the choice of the 
flavour basis becomes a delicate issue. In fact, it turns out that Boltzmann equations 
written in different bases cannot be transformed into one another by a flavour rotation 
i —>• Ui, and give different numerical results for the baryon asymmetry. This is illustrated 
in Fig. [21 which compares the B — L asymmetry computed in the flavour-covariant forma- 


7 














lism using the density matrix with the naive computation involving Boltzmann equations 
for the flavour asymmetries A i a , written either in the charged lepton or in the neutrino 
mass eigenstate basis. Since physics should not depend on a choice of basis, we are led 
to conclude that the density matrix formalism is the only valid approach when charged 
lepton Yukawa interactions are out of equilibrium^. 

In the intermediate temperature range where the tan Yukawa coupling is in equilibrium 
but the muon and electron Yukawa couplings are not, the quantum correlations between 
the tau and the other flavours are destroyed. In practice, this means that the off-diagonal 
entries in the third line and the third column of the density matrix are driven to zero by 
the fast tau Yukawa interactions. Thus, the relevant dynamical variables in this regime are 
the 2x2 density matrix (A^) Q « describing the asymmetries stored in the lepton doublets 
£ e , in and their quantum correlations, and the asymmetry A^ t stored in l T . 


3.2 The Boltzmann equation for the density matrix 

The derivation of the evolution equation for the density matrix is not straightforward. In 
this paper, we shall use the closed time path (CTP) formalism [46] . which has been used 
to obtain flavoured quantum Boltzmann equations for the standard leptogenesis scenario 
with heavy Majorana neutrinos HU [12, [T3J E, 13 |T6, 33 CEJ (for other approaches, 
also in the framework of the standard leptogenesis scenario, see e.g. Refs. 0 02 ] or 
the review ED- In this formalism, which is well adapted to describe non-equilibrium 
phenomena in quantum field theory, particle densities are replaced by Green’s functions 
defined on a closed path C in the complex time plane going from an initial instant t = 0 
to t = Too and back. When applied to leptogenesis, this formalism leads to quantum 
Boltzmann equations involving memory effects and off-shell corrections m3; in particular, 
the CP asymmetries are functions of time and their values depend on the history of the 
system. Such effects can be important for resonant leptogenesis and soft leptogenesis m 
113, but they are not expected to play a significant role in the scenario studied in this 
paper, which does not involve degenerate states. We shall therefore ignore them and 
make several simplifying assumptions in order to obtain a (classical) flavour-covariant 
Boltzmann equation for the density matrix. In particular, we shall ignore plasma/thermal 
effects and apply the CTP formalism to quantum field theory at zero temperature. This 
procedure is going to provide us with an equation of the form: 

sHz = (|| - x ) i d £ <*p - ’ ( 3 - 2 ) 

where the right-hand side contains a source term proportional to the flavour-covariant CP- 
asymmetry matrix £ a p and washout terms associated with inverse triplet and antitriplet 
decays (W^), 2 —>• 2 scatterings involving leptons and Higgs bosons (J), 2 lepton-2 

‘In the limit where scatterings involving lepton and Higgs doublets can also be neglected, there is however 
a privileged basis in which scalar triplet leptogenesis can be described in terms of the three “diagonal” flavour 
asymmetries (Ae) aa = Ag a . Indeed, in the basis where the triplet couplings to leptons f a) 3 are flavour diagonal, 
the evolution of the diagonal entries of the density matrix becomes independent of its off-diagonal entries, as 
can be checked using Eqs. (13.21) . (13.481) . (|3.53l) and (13.541) . This “3-flavour approximation”, which is valid in the 
temperature regime where all charged lepton Yukawa interactions are out of equilibrium, is the analog of the 
single flavour approximation in the leptogenesis scenario with right-handed neutrinos. 










lepton scatterings (W^g) and 2 —>■ 2 scatterings involving leptons and triplets (W^). By 
construction, Eq. (13.21) is covariant under flavour rotations £ —>• U£, which means that 
each of the matrices £, W D , W iH , W 4 ^, transforms in the same way as the density 
matrix A^, namely as M. —>• U*M.U T . 

Before proceeding with the derivation of Eq. (13.21) . let us specify our notations. For any 
species X = £ a , H , A, we define the comoving number density Y\ = nx/s, where s is the 
entropy density, and the asymmetry stored in X by Ax = Yx — Yx = ( n .Y — n x )/ s • We 
also define the comoving number density of triplets and antitriplets Ea = (n a + n^)/s. 
The evolution of these quantities as a function of z = M&/T is governed by a set of 
Boltzmann equations. Finally, a superscript eq denotes an equilibrium density. 


3.2.1 Derivation of the Boltzmann equation in the CTP formalism 

The relevant degrees of freedom in scalar triplet leptogenesis are lepton doublets, Higgs 
doublets and scalar triplets. We will therefore need the following Green’s functions, core¬ 
sponding to all possible orderings of the fields along the closed time path C (for a review 
of the CTP formalism, see Ref. [35]): 


G ap(. x iy) = -^o.{x)Ip{y)) , (3.3) 

G afj( x ’y) = *(^(y)4(z)), (3.4) 

Gipix, y) = 0( x ° - y 0 )G^(x, y) + 8(y° - x°)G<p(x, y) , (3.5) 

y) = 0(y° - x°)G>p(x, y) + 0{x° - y°)G<p(x, y) , (3.6) 

where a, /3 are lepton flavour indices, and the f’s refer to left-handed lepton doublets, 
whereas for a scalar field 4>(x) (representing a Higgs doublet or a scalar triplet): 

G ^{x,y) = -i(0(x)0 f (j/)), (3.7) 

G ^i x ,y) = -i(0 f (y)0(x)), (3.8) 

Gj(x, y) = 0{x° - y°)G £ (x, y) + 9(y° - x°)G^{x, y ), (3.9) 

G f 0 (x,y) = 6»(y 0 -x°)G>(x,y) + 6l(x 0 -y 0 )G<(x,y). (3.10) 

The brackets (...) mean that we take the average over all available states of the system. 
One can write these Green’s functions as a single 2x2 matrix: 

■>(-(? S). 

where the plus sign refers to bosons and the minus sign to fermions. This matrix satisfies 
the following Schwinger-Dyson equation: 


G(x, y) = G°(x, y) + [ d A w\ [ d i w 2 G°(x, wi)Y,(wi, W 2 )G(w 2 , y ), ( 3 . 12 ) 

Jc Jc 


where S is a 2 x 2 matrix containing the self-energy functions S > , E < , E* and E 4 , defined 
in an analogous way to the Green’s functions G > , G < , G* and Gh 


E B 

F 




( 3 . 13 ) 
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and G° is the free 2-point correlation function. The Schwinger-Dyson equation can also 
be written as 


G(x,y) = G°(x, y) + / d A w i / d A w 2 G(x, w 2 )G°(w 2 , y). 


(3.14) 


In the fermionic case, we note for later use that acting on Eqs. (13.121) and (13.141) with the 
operators i$ x and i0 y , respectively, gives the following equations of motion: 

WxG a p{x,y) = 8^{x-y)8 a pl + V' [ d A w E a7 (x, w)G y p(w, y) , (3.15) 

7 Jc 

G a p(x,y)i0 y = -5^\x-y)5 a pl + E / d 4 w G ai (x, y) , (3.16) 

7 JC 

where we have restored the flavour indices and used the fact that the free Green’s function 
for massless fermions satisfies i$ x & a p{x, y) = (x — y) 5 a p /, with I the identity matrix 
in both spinor and CTP spaces. 

Writing the free Dirac field as 


ip(x) = 


d 3 p 


(27t) 3 a /2 


E ( u (p» «)&& s ) e ~ w ' x + «(p> s ) d] & s V p ' x ) > (3-1-7) 


where u ip = y/p 2 + m 2 , we define the phase-space distribution functions of lepton and 
antilepton doublets p a / 3 {p) and p a p{p) as matrices in flavour space by 

( b Up) b p(p')) = (2n) 3 8 ( - 3) (p-p , )p a p(p) , (3.18) 

(d}p(p)d a {p')) = (2 tt) 3 5 (3) {p- p')p a p{p). (3.19) 

The reversed order of the flavour indices a and f3 in the definition of p a p{p) ensures that 
the distribution functions of lepton and antilepton doublets transform in the same way 
under a rotation in flavour space, i —>• Ut : 

p -> U*pU T , (3.20) 

p -> U*pU T . (3.21) 


Similarly, the phase-space distribution functions of a charged scalar 0 and of its antipar¬ 
ticle are defined by 

( a \{p)H(p')) = (27r) 3 <5 (3) (p-p')MP) > (3-22) 

( b l(p)Mp')) = (2vr) 3 -5 (3) (p-p')^(p) > (3-23) 

where a^(p) and 6^(p) are the annihilation and creation operators appearing in the defi¬ 
nition of the free charged scalar field: 

«» - / + ■ < 3 - 24) 
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With these definitions, the Green’s functions for left-handed lepton doublets can be 
written as (neglecting lepton masses) 


<s>. 

£ V 

Xa 

TT 

II 


d3 p f 

(:27t) 3 2o ;p * t 

[Sap - ppa(p)\ e 

- ip - {x -y) + Pp a (p ) e i p< x -y)}p R , 

(3.25) 

iG ap( X ’V) = 

-Pl 

f d 3 p 

1 ( 27t) 3 2w p y 

>{pp».(p)e- ip - (x - 

~ V) + [$ap ~ Ppa(p)\ 

e ip-(x-y) | p R 








(3.26) 

whereas for scalars one obtains 





( x - 

<y) = 

r d 3 p 

J (27r) 3 2wp- 

{[1 + Mji)]e-* 

ix ~ y) + Pp(p)e W<X ~ 

- y) } , 

(3.27) 

iGf k (x. 

<y) = 

f d 3 p 

J (2tt) 3 2uj£ 

{po{p)z- iHx ~ y) 

+ [1 +pm^ iHx ~ 

-y) | 

(3.28) 


Strictly speaking, these expressions are valid for free Green’s functions only, but they will 
be sufficient for our purpose. Notice that, under charge conjugation, 


CG>pC-\x,y) = CG<l{y,x)C~ l ee ~G< aR (y,x), (3.29) 

CG>C~\x,y) = G<(y,x), (3.30) 


where C, not to be confused with the charge conjugation operator C, is the charge conju¬ 
gation matrix defined by C'yJC '” 1 = — 7 ^, and the quantity G^ aR in Eq. (|3.29 j) is defined 
as Gi in Eq. (13.261) with Pj J interchanged with P R . 

The comoving number densities of scalars and left-handed leptons, and of their an¬ 
tiparticles, are given by 




np 

s 




Y a p = 


_ 


Y a p = 


^a/3 


d 3 p 

(27r) 3 

d 3 p 

(27r) 3 


Pap ip) , 


Pap(p) , 


(3.31) 

(3.32) 


where Y a p and Y a p, similarly to the phase-space distribution functions p a p(p) and p a p{p), 
are matrices in flavour space, and g<f, = 2 (resp. g$ = 3) for Higgs bosons (resp. scalar- 
triplets). One also defines the matrix of lepton asymmetries (hereafter called “density 
matrix”): 


(A ,) a/3 ee = Y af) - Y afi . (3.33) 

The diagonal entries of (A() a p correspond to the flavour asymmetries stored in lepton 
doublets, while the off-diagonal entries encode the quantum correlations between the 
different flavour asymmetries. From the definition (13.331) and from Eqs. (13.201) and (13.211) 
one can see that Ai transforms under a rotation in flavour space i —>• Ut as 


-»• U*A e U T . 


(3.34) 
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One can show that An Q( g is the zeroth component of the current : (or 

more precisely of its average (,7^)): 

(3.35) 


(3.36) 


A n aP = {J%) = (: v a ip :). 

We will use this fact to derive an evolution equation for A n a g. Noticing that 


( d v J *p) = ~ tr {iVx + i$y)G} a (x,y) 


iy = x 


where the trace is taken over spinorial and SU (2 )l indices, and using the Schwinger-Dyson 
equations (13.151) and (13.161) to express the right-hand side of Eq. (13.361) in terms of the 
self-energy functions E > and E < , one obtains 

( d V J ap) = ~ [ d4 ™ tr ^( x ^ w ) G ya( W ^ X ) 

J c, 

-G^ 7 (x, w^^w, x) + G^ 7 (x, w)Y,> a (w, x) 


= — d 3 w / dt w tr 


Y,> 7 (x,w)G< a (w,x) - ^ 7 (x,w)G> a (w,x) 


-G| 7 (x, w)T,< a (w, x ) + G| 7 (x, w)Y,> a (w, x) 


(3.37) 


Since we consider a homogeneous and isotropic medium, the divergence of the current 
reduces to dAn a g/dt. Finally, we incorporate the expansion of the universe by making 
the following replacement in the above equation: 




dAn a p 


+ 3HAn a i3 = sHz 


d{Af>') a ig 


dt p dz ' 

We thus obtain the quantum Boltzmann equation for the density matrix (Af_) a g: 

d{A^) a g 


(3.38) 


sHz- 


dt w tr 


^( x ,w)G 7 a (w,x) - ^ 7 (x,w)G> a (w,x) 


&e) a p f f 

-Gp 7 (x,w)T,< a (w,x) + Gp 7 (x,w)^ 7a {w,x) 


(3.39) 


This equation is both quantum and flavour-covariant, but we are only interested in flavour 
effects. We shall therefore take the classical limit by keeping only the contributions to E 
involving an on-shell intermediate part. This procedure will provide us with a flavour- 
covariant Boltzmann equation of the fornil (13.21) . 

8 Strictly speaking, the Boltzmann equation for the 3x3 density matrix (Ae) a p is valid only in the regime 

where the quantum correlations between the different lepton flavours are not affected by the charged lepton 
Yukawa interactions, i.e. at T > 10 8 * * * 12 GeV. When this condition is not satisfied, one must either add a term 

accounting for the effects of the Yukawa-induced processes on the right-hand side of Eq. <& or impose that 

the appropriate entries of the density matrix A^ vanish (see discussion later in this section and Section [5] for the 

relevant Boltzmann equations). Furthermore, the effect of spectator processes such as sphalerons and Yukawa 
interactions, which impose relations among the various particle asymmetries in the plasma, is not included at 
this stage (it will be discussed in Section |4|. 
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ra.r~ x ('r 

w “ v~ > a ) 

Figure 3: 1-loop contribution to the lepton doublet self-energy Y,p a (x,y). 

3.2.2 Washout terms: decays and inverse decays 

Let us first compute the flavour-covariant washout term associated with triplet/anti- 
triplet decays and inverse decays. In the CTP formalism, this term arises from the 1-loop 
contribution to the lepton doublet self-energy shown in Fig. [3j For the first term of the 
integrand on the right-hand side of Eq. (13.391) . this gives 


tr 


^(x,w)G< a (w,x) 


= (-*//Jp)(-*/V) tr 

= 3 ./ (2 J (2ir) 3 2uj£ / (27T) 3 2o; r 2(U) 

x [pa{p) e-^-^ + [1 + p A (p)} 
x { Ppa {k) e ~ ik -( w ~ x ) + [ Spa _ p pa (k)} e ik < w ~ x )} 
x I p n Jl) e~ u < w -*) + [5 ai - p ai {f)} e il -( w ~ x ) J , (3.40) 


where the factor 3 comes from the trace over SU{2)i indices. The integration over the 
spatial coordinates of w gives a momentum-conserving delta function. In the integral 
over t w , we take the classical limit by making the usual assumption that t is much smaller 
than the relaxation time of particle distributions, which can therefore be factorized out 
of the integral, but much larger than the typical duration of a collision, so that the time 
integral can be extended to infinity. This amounts to replacing the integral of oscillating 
exponentials by energy-conserving delta functions, plus terms proportional to the principal 
values of l/(up±uj^±u}^), where p, k and l are the momenta of the scalar triplet and of the 
two leptons, respectively. The latter terms, however, can be neglected because they arise 
at second order in the CP asymmetry. Indeed, all terms involving a principal value are 
proportional to r) p<J (k ) + fj pa (k), where p pa (k) (resp. fj pa [k)) parametrizes the departure 
of the phase-space density p ptT (k ) (resp. p pa [k)) from its equilibrium value: 


P P a(k) = p e e q (k) 


3ocr T Vocrik) 


= n e 3(h\ 




(3.41) 


Since the unbalance between the lepton and antilepton densities is generated by the asym¬ 
metries in triplet decays, which are small numbers of order e, while the lepton and an- 
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tileptonr 
one haq£| 


'opulations are maintained close to equilibrium by fast electroweak interactions, 


ripaik) - rja P (k) = 0(e), (3.42) 

r]pa(k) + fjcrp(k) = £>(e 2 ). (3.43) 


Therefore, the terms proportional to rj pa (k) + rj pa (k) on the right-hand side of the Boltz¬ 
mann equation (13.3911 are of order e 2 and can safely be neglected. By doing so one keeps 
only the terms that are kinematically allowed for on-shell particles. Dropping Bose en¬ 
hancement and Pauli blocking factors, we obtain 


IH 


dt w tr 


E%y(,x,w)G< a (w, x) 


~ / Wjkr f / (2^ / p -k-i) 

X ^ P p<J (k ) Pct'y (1) 4“ P A (p) $ P (T ^cry J' • 


(3.44) 


Proceeding in the same way to compute the other contributions from Fig. [3] to the right- 
hand side of the Boltwmann equation (13.3911 . and introducing the space-time density of 
triplet and antitriplet decays: 


7 D = 


f d3 P 

f d 3 k 

f dH 

1 (27t) 3 2 up] 

(2vr) 3 2o;^ J 

(27r) 3 2w £ - 


3 (A 2 + \ 2 H ) (k.l) 


x ( 2ir) 4 5 {A \p -k-l ) {p c £(p) + p e £{p)} , 
we obtain the washout term associated with decays and inverse decays: 

2Bi (// f )c^ +^^(YfY T f^ + fY T fY-Y^Y) 


= 


A? L 


/ a/3 


(3.45) 


ID, (3.46) 


where we remind the reader that A a = (n* — n^J/s and £a = (ba + )/ s - One can 

linearize this expression using again the fact that flavour-blind gauge interactions keep 
the lepton densities close to their equilibrium values: 


Y a fi Y a g — (A^)q.^ , 

Y a p T Pq:/ 3 — 2y; q [5 a0 + 0(e 2 )] , 


which hnally gives 


wS = 


2 B e 
A 2 


(f f )o;/5 ^-,eq 


(2/Af/t + //tA, + A,//t) 


a/3 


7 D • 


(3.47) 


(3.48) 


Note that the washout term (I3.48P contains a piece proportional to Aa, which is due 
to decays. This term appears on the right-hand side of the Boltzmann equation for 

9 Eqs. (13.421) and (13.431) generalize the relations between equilibrium number densities ne a — ni a ~ 4 Pe a T 2 
and ne a +nj> = 2 n e e q (pe = 0) [l + O (pe a /T) 2 ] , where pe a is the chemical potential of the lepton flavour a and 
p t JT = 0(e). 
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the 3x3 density matrix (A £ ) a / 3 , Eq. (13.211 . One can easily check that it transforms as 
W D —» U*W D U T under flavour rotations £ —>• U£, as required by flavour covariance. 
Below T = 10 12 GeV, however, the tau Yukawa coupling is in equilibrium and drives the 
(e, r), (/U, r), (r, e) and (r, /r) entries of the density matrix to zero. The relevant dynamical 
variables in this regime (or more precisely in the temperature range 10 9 GeV < T < 
10 12 GeV, before the muon Yukawa coupling enters equilibrium) are A £t , the asymmetry 
stored in £ r , and a 2 x 2 matrix A 9 describing the flavour asymmetries stored in the lepton 
doublets £ e ,£n and their quantum correlations. The corresponding washout terms VV^ 
(where a and f3 label any two orthogonal directions in the (£ e , £^) flavour subspace) and 
VV ^ 2 are simply obtained by setting (A £ ) aT = (A|) TQ = A £r S aT in Eq. (|3.48p . yielding 


W d b = — 

y v a/3 ^2 


{fP)a0j£ 


+ 4 yeq ( 2 /(A 9 ) T / 1 ' + fp A 9 + A eff^ a p + 2 y e q f a Tp 


■A 


and 


W D = 

T Xj 


(fP) rr^+ 1 


TT 2 Y eq 


A p T P )tt + ((/ P)tt + l/rr| 2 ) A £j 


ID , (3.49) 


ID- (3.50) 


The resulting Boltzmann equation for (A 9 ) a( g is covariant under flavour rotations in the 
(£ e , £p) subspace. Finally, below T = 10 9 GeV, the muon Yukawa coupling enters equilib¬ 
rium and drives the (e, /r) and (fi, e) entries of the density matrix to zero. The Boltzmann 
equation (13.21) then reduces to three equations for the flavour asymmetries A £a = (A £ ) aa 
('a = e, n, r), with a washout term W® given by 


< = Y i uf 

A l 


Aa A^ a + A. 
+ 2Y/ q 


b 


7 D ■ 


(3.51) 


3.2.3 Other washout terms: 2 —> 2 scatterings 

The other washout terms in the Boltzmann equation (13.21) correspond to 2 —» 2 scatterings 
and arise from 2-loop contributions to the lepton doublet self-energy E^ a . The term VV^g 
accounts for the washout of the flavour asymmetries by the A L = 2 scatterings £^£g <->• HH 
and £^H £$H and is given by 


<8 = 


2 7 


+ 


(2/Aj/t + / pAe + A £ f P) a p A Hft ^ 

4 Y/ q + v e q [ff) 

(2fAfK f + fK f A £ -hA £ fK f ) a p 


1m 


1 


+ 


+ 


Re [tr(/«t)] 

1 

Re [tr(//ct)] 

1 


4Y/ q 

(2kAJ p + KpAg + AgKp) a p Ah 


+ yeq (f^)aP 


4Y eq 


+ yeqPP)a/3 


ItH 


7 £H 


A l 


(2kAJV + KK tA £ + A £ kk)) Ah 
4Y eq ^ 


7£H 


(3.52) 
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where A h = (n# — n#)/s is the asymmetry stored in the Higgs doublet, X K = ^/tr(«uct) 
and 7 ^-, 7 ^ and 7 ^ are the contributions (summed over lepton flavours) of different 
self-energy diagrams to the space-time density of A L = 2 scatterings ')ih- Namely, one 
has 7 w = 7 iH + 27 f H + 'yJjj, in which 7 ^ is the contribution of scalar triplet exchange, 
7 Jjj is the contribution of the D = 5 operator (12.91) responsible for m%, and is the 
interference term. Explicit expressions for the reduced cross-sections that are needed to 
compute numerically these reaction densities’ll (as well as the ones appearing in Eqs. (13.531) 
and (13.541) below) can be found in Appendix [Aj 

The washout term due to the 2 lepton-2 lepton scatterings l p l a and l^t p -H- igt a 

reads 


yy 4 ( = _ 

r y a /3 ^4 


(2/A7t + //tA f + A f //t) a)! tr(A ,//*) ,,,„ 

z^Y' ec \. v' ec i \J j )ot.fi 


yp 


lAl- 


(3.53) 


whereas the washout term associated with the lepton-triplet scatterings Ay A -H- £ 5 A, 
Ay A -H- A$A and 1^1$ -H- AA is 


W^ = 


tr(//t//t) L 2 E/ q 


(//t//tA f - 2//tAJ/t + A,//t//t) 


a/3 


7£A • 


(3.54) 


This completes the derivation of the washout terms appearing on the right-hand side of 
Eq. (13.21) . The expressions of W^jg, and relevant in the temperature ranges 

10 9 GeV < T < 10 12 GeV and T < 10 9 GeV can be obtained from Eqs. (13.521) . (13.531) 
and (|3.54l) following the same procedure as the one used for , see the discussion above 
Eqs. (|3.49l) and (13.511) . Note that and vanish in the single flavour approximation 
where a and (3 take a single value; this is consistent with the fact that the 2 lepton -2 
lepton and lepton-triplet scatterings affect leptogenesis only when lepton flavour effects 
are taken into account. 


3.2.4 Source term 

Finally, let us derive the source term of the Boltzmann equation (13.21) . In the CTP 
formalism, the 2-loop self-energy diagrams shown in Fig. [4] when inserted in the right- 
hand side of Eq. (13.391) , give rise to the flavour-covariant source term 


Sa/3 — 



7 D £a/3 , 


(3.55) 


where £ a p is a Hermitian matrix transforming as £ —>• U*£U T under a flavour rotation 
t —>• Ui. The explicit form of £ a p can be obtained without going through the full CTP 
computation by noticing that, in the temperature regime T < 10 9 GeV where the quantum 
correlations among lepton flavours are destroyed by fast interactions induced by the muon 

10 For some reactions the contribution of real intermediate states must be subtracted in order to avoid double 
counting in the Boltzmann equations. Although this subtraction procedure is embedded in the CTP formalism 
iself [T3] [15] , in practice we first take the classical limit in the quantum Boltzmann equation in order to obtain 
Eqs. (13.521) . (13.53|) and (13.541) . and then insert the subtracted reaction densities computed in AppendixlAlinto 
these expressions. 
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Figure 4: 2-loop contributions to the lepton doublet self-energy E p a giving rise to the CP 
asymmetry £ a p. 


and tau Yukawa couplings, the source term for the flavour asymmetry Ai a = (A g) aa 
(a = e, /j,, t) reads 


£rvrv — 



(3.56) 


where e Q7 , not to be confused with £ ai , is the CP asymmetry in the decay A —>• £ a £ry 
defined in Eq. (12.1311 . Comparing Eq. (|3.55l) with Eq. (|3.56[) and using Eq. (12.141) . one 
obtains the diagonal entries of £\ 


f — 


E 




1 M a 
47r v 2 


Im 


V B^Bh 


(mu m\ 


A J aa 


m a 


(3.57) 


It is then easy to derive the full flavour-covariant CP-asymmetry matrix £ by reading its 
dependence on the couplings f a p and n a p from Fig. |4l and by using its hermiticity and 
covariance properties together with Eq. (|3.57D : 

£ a0 = 4-4^ SB ,. (3.58) 

om v A rriA 

The first equality in Eq. (13.5711 implies that the trace of this matrix is equal to the total 
CP asymmetry in triplet decays: 


tr£ = ^2^a/3 = e A • (3.59) 

a,(3 

This completes the derivation of the right-hand side of the flavour-covariant Boltzmann 
equation (13.21) . 


3.3 Full set of Boltzmann equations 

Finally, Eq. (13.21) must be supplemented with the Boltzmann equations for Ea and Aa 
( an equation for A h is not needed since, as we are going to see in Section |4j A h can be 
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expressed as a function of the other asymmetries): 


sHz it=-a t ■ - 2 ((it) ■ ') 74 ’ (3 - 6o) 

sHz ^hl = (|| - l) Id S a p - W& - W- - Wja - , (3.61) 

~ (tr(W D ) - Wg) , (3.62) 

where the first and second terms in Eq. (13.601) come from (anti)triplet decays and triplet- 
antitriplet annihilations, respectively, and Wg in Eq. (13.621) is the washout term due to 
the inverse decays HH —>• A and HH —> A: 

Wg = 2B H (^-^y D . (3,63) 

Using Eqs. (|3.48l) and (13.631) . the Boltzmann equation for A a can be rewritten as 


sHz 


dA& 

dz 


( Aa , „ tr(//tA<) 

Vs e A q A 2 v; q 


R A H \ 

Bh y^ ) 1d 


(3.64) 


Eqs. (13.601) and (13.621) take the same form as in the single flavour case (261 , except for the 
term associated with the decays A -H- £*,£$, A -H- £^£§ in the Boltzmann equation for Aa, 
which involves a trace with the density matrix rather than the total lepton asymmetry. 


4 Spectator processes 

Various Standard Model (SM) reactions (strong and electroweak sphalerons, Yukawa cou¬ 
plings) that contribute indirectly to the washout of the lepton flavour asymmetries have 
to be included in the computation. Indeed, while these processes cannot create or erase a 
lepton asymmetry by themselves (for this reason, they are usually dubbed spectator pro¬ 
cesses), they modify the densities of the species on which the washout rates depend, thus 
affecting the final baryon asymmetry. For instance, Ac t can be turned into A h and A Tr 
by y r -induced interactions, and into A 93 by electroweak sphalerons (where, as previously 
defined, we denote the asymmetry stored in the species X by Ax = (nx — ^x)/ s )- 

In this paper, we follow the standard approach by assuming that each spectator process 
is either negligible or strongly in equilibrium, in which case it imposes a relation among 
the chemical potentials (hence among the asymmetries) of the species involved. The set 
of relations valid in a given temperature range is then used to express the asymmetries 
appearing in the Boltzmann equations in terms of asymmetries that are conserved by all 
active SM processes. This was first done in the standard leptogenesis scenario with heavy 
Major ana neutrinos in Refs. M, and in scalar triplet leptogenesis in Ref. m- In the 
following, we generalize the analysis of Ref. [31] to the density matrix formalism. 

The asymmetries appearing in the Boltzmann equations (|3.60l) - (13.62|) are Aa, A# and 
the density matrix (A^) a ^ (which reduces to A^ r and the 2x2 matrix (A q ) a ^ in the range 
10 9 GeV <T < 10 12 GeV, and to the flavour asymmetries A^ e/iT below T = 10 9 GeV). 
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Out of these, only Aa is conserved by all SM interactions, while (A n) a p is affected by 
electroweak sphalerons and charged lepton Yukawa couplings. It is therefore convenient 
to introduce the following flavour-covariant matrix: 

A Q /3 = g 5 a p — (A e ) a p , (4.1) 

which is conserved by all SM processes except the charged lepton Yukawa interaction^], 
and the asymmetries: 

= &B/ 3 -L a = g - A £a - A eRa (a = e, p, r), (4.2) 

which are preserved by all SM interactions. Above T > 10 12 GeV, where electroweak 
sphalerons and charged lepton Yukawa couplings are out of equilibrium, A a p is conserved 
and reduces to —(Ag) a p. Below T = 10 9 GeV, the tau Yukawa coupling is in equilibrium 
and Ai t is partially converted into A TR . The conserved asymmetries in the temperature 
range 10 9 GeV < T < 10 12 GeV are therefore 

A ° a/3 = -A b 8 a p - (A 9 )^ , A r and Aa , (4.3) 

where a and f3 label any two orthogonal directions in the (£ e , l R ) flavour subspace. Finally, 
below T = 10 9 GeV, the muon Yukawa coupling is also in equilibrium and A 9 ^ must be 
replaced by A e and A^. 

Using chemical equilibrium, one can express the Higgs and lepton flavour asymmetries 
as functions of the conserved asymmetries in each temperature range. Thus, above T = 

10 12 GeV: 

{Ai)a/3 = ~A a /3, = g H (A pa , Aa) , (4.4) 

while between T = 10 9 GeV and T = 10 12 GeV: 

(A °) q/ 3 = g^p(A° pa , A t , Aa) , A, r = g e T (A° pa , A r , A a ) , A H = g H (A° pa , A r , A a ) , (4.5) 

and below T = 10 9 GeV: 

A £q = gi( A p ,A a ), Ah = g H { A p ,A a ), (4.6) 

where the functions g e and g H and their arguments depend on the temperature. 

The relations among chemical potentials associated with the various spectator pro- 

11 The diagonal entries of A a p reduce to iA s - A £q , which is conserved by electroweak sphalerons. Its 
off-diagonal entries are conserved by electroweak sphalerons as well since there is no flavour off-diagonal lepton 
anomaly (d p (j£) a p = 0 for a ^ /3). 
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cesses are (we assume gauge interactions to be in equilibrium) 


Udi) 

- Y (Wa+MeJ 

I + Y 6/i A — 0 (Ay — 0) 

(4.7) 

i=l,2,3 

a=e,fi,r 



(^1% + AGj + l^di) 

i=l,2,3 

= 0 

(A b = 0 ) 

(4.8) 

Y ( 2 ^« “ - Md*) 

*=1,2,3 

= 0 

(QCD sphalerons) 

(4.9) 

Y + Y ^ 

= 0 

(EW sphalerons) 

(4.10) 

*=1,2,3 a=e,/j.,T 




~ Hut Y fJ'H = 0 


(up-type quark Yukawa) 

(4.11) 

Hqi - Udi ~ HH = 0 


(down-type quark Yukawa) 

(4.12) 

MG - Me a - HH = 0 


(charged lepton Yukawa) 

(4.13) 


in which we have also included the constraint due to the hypercharge neutrality of the 
universe, as well as Ag = 0 , which holds as long as the electroweak sphalerons remain out 
of equilibrium. These constraints on chemical potentials can be translated into relations 
among particle asymmetries using the formulae, valid at leading order in n/T\ 

A t = | w T 2 , A, = g (4,14) 

for a bosonic or fermionic species, respectively, with g & (gf) the number of internal degrees 
of freedom. 

Taking into account the processes that are in equilibrium in a given temperature 
rang^l, one arrives at the following expressions for the functions g e and g H : 

(i) T > 10 15 GeV. In this temperature range all spectator processes are out of equilib¬ 
rium (but equality of the chemical potentials within the same gauge multiplet is still 
assumed), so one simply has, from Ay = Ab = 0: 

(A t) a p = —A a/ 3, (4.15) 

A h = —tr(A a/ g) - 2 A a . (4.16) 

(%%) 10 13 GeV < T < 10 15 GeV. Only the top quark Yukawa coupling and gauge inter¬ 
actions are in equilibrium. Asymmetries in and t B can now develop, but baryon 
number is still conserved. The constraints (14.711 . (14.8p and (14.111) give 


(A^) Q( 3 = —Aa/3, (4.17) 

A// = —-tr(A Qi a) —-A a • (4-18) 

12 Using the reaction rates available in the literature [151 150115T1 152j . one finds that the top Yukawa coupling 
comes into equilibrium around T = 10 15 GeV, the QCD sphalerons around T = 10 13 GeV, the bottom and 
tau Yukawa couplings as well as the electroweak sphalerons a little beneath T = 10 12 GeV, the charm Yukawa 
coupling around T = 10 11 GeV, the strange and muon Yukawa couplings around T = 10 9 GeV, and the electron 
Yukawa coupling around T = 10 5 GeV. 
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(in) 10 12 GeV < T < 10 13 GeV. In this range, the QCD sphalerons are also in equilib¬ 
rium. Asymmetries in all quark species can develop, but neither baryon number nor 
quark flavour are violated, hence A^ = A b 2 = Ar 3 = 0. Moreover, the first two 
generations of quarks and bn are created only through QCD sphalerons, implying 
A qi = A q2 = —2A U = —2Ad = — 2A S = — 2A C = —2A&. This together with the 
constraints (14.71) (14.911 and (14.111) yields 

{&t)aP = —AaP, (4.19) 

14 28 

Ah = _ 23 tr ( A "' S ) ~ 23 Aa ' ( 4 - 20 ) 

(iv) 10 9 GeV < T < 10 12 GeV. In addition to the previous processes, the electroweak 
sphalerons and the bottom, tau and charm Yukawa couplings are also in equilibrium 
(we neglect the difference between the charm and bottom quark equilibrium tem¬ 
peratures). Quantum correlation between l T and the other lepton doublet flavours 
disappear, and the relevant asymmetries in the lepton sector are A 9 ^ and A r . Due 
to electroweak sphalerons, equal asymmetries develop in each quark generation: 
A b 1 = A b 2 = ^b 3 (/ 0). Moreover, since the right-handed quarks ur, cIr and 
sr are created only through QCD sphalerons, one has A u = A^ = A s . Together 
with Eqs. (14.71) and (14.91) (j4. 131) . this constraint gives 

(A ?) a/3 - tr(A° /3 ) + — A t + 5 a p - A% , 

A, = —tr(A 9 ) - — A r - — A a , 

589 V af)1 589 589 

. 164 , a0 , 224 A 344 A 

A h =-tr(A„o)-A r - A a . 

589 v afiJ 589 589 


(4.21) 

(4.22) 

(4.23) 


(v) 10 5 GeV < T < 10 9 GeV. In this range, the second generation Yukawa couplings 
are also in equilibrium. All three lepton flavours are now distinguishable and the 
relevant asymmetries in the lepton sector are the A a . The QCD sphalerons lead to 
A u = Ad- Together with this constraint, Eqs. (14.71) and (14.91) (14.131) yield 


A 4 

A4 

A h 


151 . 20 A 20 A 4 A 

-A e H - - Am -\- -A r +- Aa 5 

179 179 M 179 179 

25 A 344 A 14 A 11 

358 537 M 537 179 

25 A 14 344 11 

-A e H- A a - A t -A a , 

358 537 M 537 T 179 ’ 

37 A 52 A 52 A 82 A 

-Ae-A u -A r - A a . 

179 179 M 179 179 


(4.24) 

(4.25) 

(4.26) 

(4.27) 


(vi) T < 10 5 GeV. All spectator processes are now in equilibrium. Eqs. (14.71) and (I4.9D - 
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(I4.13P then give 


a 4 

a 4 


Ah 


442 A 32 A 32 A 3 A 

-Ag -A n ~\~ -A T — —Aa 5 

711 711 M 711 79 

32 A 442 A 32 A 3 A 

-Ag —-A n “h-A r — —Aa i 

711 711 M 711 79 

32 A 32 . 442 A 3 A 

-A e H-A„-A r - A a , 

711 711 M 711 79 

16 A 16 A 16 A 26 A 

-A e -A u -A r - Aa • 

79 79 M 79 79 


(4.28) 

(4.29) 

(4.30) 

(4.31) 


5 Boltzmann equations 

Using the results of Sections [3] and [dj one can write the full set of Boltzmann equations 
for scalar triplet leptogenesis including all relevant spectator processes in a covariant way. 

The dynamical variables are Ea = (ka + n^)/s, Aa = { n A ~ n /\)/ s an d, depending on 
the temperature regime, the 3x3 density matrix A a p, the 2x2 density matrix A^ and 
the asymmetry A T , or the three asymmetries A e , A^ and A r . 

At temperatures higher than 10 12 GeV, the system of Boltzmann equations i J^I 

= - 1 ) 7 - 4 - (5J) 
sHz ^ “ ^jlD £a/3 + + wg, + Wij , (5.2) 

sHz ^r = ~\ ( tr ( wD ) - w h ) > (5-3) 

with the flavour-covariant CP-asymmetry matrix £ a p given by Eq. (|3.58l) . and the washout 
terms W^, W^, and given by Eqs. (13.4811 . (13.521) . (I3.53H . ([3. 541) and 

(|3.63j) . respectively. In the expressions for the washout terms, the asymmetries (A^) a p and 
A h should be substituted for —A a p and g H (A pcr , Aa), respectively, where the function 
g H depends on the temperature and can be found in Section 01 Cases (i) to (in). Once 
this is done, the system of equations (15.1[) - (15.3D has a closed form and can be solved 
numerically. 

When the temperature drops below 10 12 GeV, the system of Boltzmann equations 

13 One could write the Boltzmann equations in terms of the density matrix (Ae) a p rather than A a p, since 
(Ae) a p is conserved by all SM processes above T = 10 12 GeV. However, it is more convenient to use A a p as 
a dynamical variable, as (A t) a p is no longer preserved by electroweak sphalerons when the temperature drops 
below 10 12 GeV. 
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becomes 


sHz 


<E a 


sHz- 


dz 

dA<! 


a/3 


sHz 


dz 

dA r 


sHz 


dz 

dA A 

dz 


(s? 2 ((s|) 1 ) X4 ’ 

“ ( _ 1 ) ^ 01/3 + + ^"Z 3 + ( a ’ ^ 

~ (|| - l) ID £rr + VV r D + Wi H + Wf + M\ 

"(tr (W D ) + W?-W°), 


(5.4) 

e,ii), (5.5) 

(5.6) 

(5.7) 


where the washout terms WY (a, /3 = e, //) and >V T " can be obtained by setting (A^) QT = 
(A|) Ta = Ai T 5 aT in the expressions of the corresponding W/Y (a,/? = e,/r,r), and the 
asymmetries (A/) ai g, Ag T and A h appearing in the washout terms should be replaced by 
the functions g e a/3 (A° prT , A r , A a ), g £ T (A° p(T , A r , A a ) and g H (A° pal A r , A a ) given in Sectional 
Case (iv). Finally, below T = 10 9 GeV: 


sHz 


<E a 

dz 


sHz 


dA a 

dz 


sHz 


dA A 

dz 



where e a = £ aa = ^ 7 e a7 , the washout terms VV a ” can be obtained by setting (A^) Q( g = 
A e a da/3 hi the expressions of the corresponding W Q '/, and the asymmetries Ap a and Ah 
appearing in the washout terms should be substituted for the functions g^(A p , A a ) and 
g H (A p , A a ) defined in Section^ Cases (v) and (vi). 

In Section 0 we will compare the results obtained with the above system of covariant 
Boltzmann equations (hereafter referred to as the full computation) to various approxi¬ 
mations in which flavour effects and/or spectator processes are neglected. We give below 
the equations corresponding to these approximations. 


Flavoured computation, no spectator processes In this approximation, all 
SM Yukawa interactions as well as the QCD and electroweak sphalerons are neglected 
(rigorously, this is legitimate only above T = 10 15 GeV, but here it is assumed for arbitrary 
temperatures). The relevant Boltzmann equations are Eqs. (15.111 (15.311 with (A^)^ and 
A h in the washout terms replaced by (Case (i) of Section 0|) 

{^i)ap = ~^a/3 , A H = - tr(A a(S ) — 2 A a . (5.11) 

Single flavour approximation, no spectator processes In this approximation, 
both flavour effects and spectator processes are neglected. To derive the relevant Boltz¬ 
mann equations, we start from Eqs. (I3.60P - (|3.62D and we implement the single flavour 
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approximation by substituting 


(A-e ) a /3 —>• A^, f a p —> \i , K a p —> A k , £ a p —>■ tr£ — 6 a • (5-12) 

This gives: 

= - (If- ‘K - 2 ((If ) 2 - (A p' 13 ) 

= (|| - 1 ) 7 ^ + W t D + W eH , (5.14) 

sHz ^r = \ ( W(D ~ ’ (5 - 15) 

where and W® are the washout terms due to inverse triplet and antitriplet decays: 
W, D = -2B, (A + ^) 7D , Wg = -2B„ (Ag - ||) 7D , (5.16) 


while W (H represents the 2 —>■ 2 scatterings involving two leptons and two Higgs bosons: 


W (H = 




^ Y ®q 1 Yff-) ' 

Then we impose the single flavour version of the relations (15.111) : 

= -A , A h = —A - 2Aa , 

where A = A b-l- This gives 


dS A (T, a . 

sff2 ^r = -(E«- 1|7c - 2 


Sa\ 2 _ 

^aJ 


1 7a 


sHz ^r = 

dz 


sHz 


c/Aa 

dz 


Sa A _ 9R / A _ A a \ 
j-ieq lj7D e A y] ec l J YD 


- 2 ^A + A + 2Aa ^1 


A / 


_A_„A 

U e A q 


A + 2 Aa^\ 

+ Bh -^eq- J 7D 


7m 


(5.17) 

(5.18) 

(5.19) 

(5.20) 

(5.21) 


This case has been studied in detail in Ref. [26] . 

We stress again that the single flavour approximation is not a limiting case of the 
flavoured computation, in the sense that it does not become exact when the charged lepton 
Yukawa couplings are sent to zero. This is why, in some regions of the parameter space, the 
single flavour approximation is actually a very bad approximation to the flavour-dependent 
computation, even though charged lepton Yukawa interactions are out of equilibrium. This 
behaviour, which contrasts with the standard leptogenesis scenario involving hierarchical 
heavy Majorana neutrinos, will be demonstrated numerically in Section [ 6 l 
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Flavour non-covariant approach, spectator processes included In this case, 
we take into account spectator processes but do not require flavour-covariance of the Boltz¬ 
mann equations. We assume instead that, in each temperature regime, a basis for the la s 
can be chosen in which scalar triplet leptogenesis is described by the evolution of “diago¬ 
nal” lepton flavour asymmetries. This is the approach that has been followed so far in the 
literature [30l I3T] (however without taking into account spectator processes in Ref. [30)). 
and it is known to be a good approximation in the standard leptogenesis scenario with 
heavy Majorana neutrinos [8], except at the transition between two temperature regimes. 
In practice, one works within the single flavour approximation above T = 10 12 GeV, and 
in a two-flavour approximation in the temperature range 10 9 GeV < T < 10 12 GeV, where 
the fast tau-Yukawa interactions destroy the coherence between the tau and the other two 
flavours. Finally, below T = 10 9 GeV, all three lepton flavours are distinguishable and 
the relevant dynamical variables are the flavour asymmetries A a , a = e,/r,r; in this case 
the Boltzmann equations are given by Eqs. (I5.8I) - (|5.10D . 

Thus, for T > 10 12 GeV, the Boltzmann equations are 


sHz 


dT, A 

dz 



sIiz ir = 

dz 


sHz 


dA A 

dz 


Sa\ 2 _ 


1 \1A, 


/£A \ „ / A A a \ 

_ y XI eq _ ^ — 2 Dt ^ ye q — ^eq ) 1/D ■ 


/ 


/A g H { A,A a ) \ 


U eq 


Y? 


JllH ■ 


Bh 9H[A ’ Aa) ^ 


V S A 


Yr 


Y? 


7 D 


(5.22) 


(5.23) 

(5.24) 


where g H ( A,A a ) = —(A + 2A a ), —2(A + 2A a )/3 and —14(A + 2A a )/23 for T > 
10 15 GeV, 10 13 GeV < T < 10 15 GeV and 10 12 GeV < T < 10 13 GeV, respectively. For 
10 9 GeV < T < 10 12 GeV we work in a 2-flavour approximation, as was done in Ref. m- 
Namely, we consider only the asymmetry stored in l T and the overall asymmetry stored 
in the l e , 1^ doublets, i.e. we perform the single-flavour approximation in the (l e , 1^) 
subspace by substituting 


where a, f3 


(A®) a p —> Ap 0 , A 9 ^ — > Aq , ~ > tr £' 0 = eo , 


''a. (3 


fa.fi ->• \J tr(//f) = /oo , far y/\fer \ 2 + \f^r \ 2 = far , 

V\ 


K a0 ->• yj tr(«/«t) = Koo , K, 

e, \i and 


ACer | 2 T | AC^t | 2 — ACqt 


£° = 


£ee £, 


efi 


'fie 


fifi 


f = 


fee fefi 
f\ie fjiji 


K = 


AC, 


qte 


e/i 


(5.25) 


(5.26) 
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We thus obtain the following system of Boltzmann equations (in which a, 7 = 0,r): 


sHz 


dz 

dA 






— 1 7-D 2 eq 1 7 A , 


S A 




(5.27) 


sHz ~JT = - + 2 B e^2 c , 


dz 
+ 2 


-iD 

y a.') 


'gi( A P, a a ) + ^( a p> a a) , a a \ 

+ ^a 7D 


2 Y/ q 


s e A q y 


( 9gA P i Aa) + 9y(A p , Aa) g H (A p , Aa) \ ,„a a , r z t , ft \ 

/ y [ 9^ e< l ' yxcq / K^ajilH ' ^arfllH ' ^a^filH) 


7 


sT/T 


+E 

7 

dA A 

dz 


2 y/ q 

'^(Ap, A a ) -4 ( a p,A a ) 

y/ q 


(c^ 7 « + 0 <a)> 


(5.28) 


Aa , D + ^(Ap, Aa) g H ( A p ,Aa)\ 

^eq + ^ 2^ C V^-^TTTea-7% 


a, 7 


2 y eq 


y h 


)7 D, 

(5.29) 




where e r = £ rT , 
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Aa , 


p / \ * \ 503 . 60 A 

ff 0 (A p ,A A ) - -—A 0 + —A r + 

g £ T (A p , A a ) = — A 0 - — A t - — A a , 

y T v p, a; 589 u 589 r 589 A 

HtK a \ 164 A 224 A 344 A 

H ' A A -'- An-— A t -— A a , 


q (A 0 , Aa) =- ^0 

y v p ’ ’ 589 589 


589 


(5.30) 


and the coefficients C^ are given by (a, 7 = 0 , r) 


/« 2 


/7L» _ JUU 

°00 — ^ 2 " > 

y^tA y^lD 

^OL'y ^ary 5 


s~iD s~iD f Or 

°o t — Wo ~ "^ 2 “ 


CL = 2 


fa 7/^1 


a r y nj a'y 


“ 7 Re[tr(/t«)] 


pO _ |/rr| 2 
U rr ^2 ’ 

■H _ Kg-y 
U Q!7 ^2 


p-, 4 £ _ p-, 4 £ _ n D 
^Or — Wd — '-'Or ) 


ZA sylA l/oo/or + /Or/r 


Wt — '- y ^n — 


r 0 


tr(//t//t) 


(5.31) 


6 Numerical study 

In this section, we investigate numerically the impact of lepton flavour effects on scalar 
triplet leptogenesis using the flavour-covariant formalism introduced previously. In order 
to assess quantitatively the relevance of lepton flavour covariance, we compare the results 
of the density matrix computation with the ones obtained in flavour non-covariant ap¬ 
proximations. We also compare the relative impacts of spectator processes and of lepton 
flavour covariance on the generated baryon asymmetry. More specifically: 

• in the temperature regime T > 10 12 GeV, where all charged lepton Yukawa couplings 
are out of equilibrium, the flavour-covariant computation involves the 3x3 density 
matrix A a g, whose evolution is governed by the Boltzmann equation (15.21) . The 
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result of this computation is compared with the one obtained in the single flavour 
approximation, for which the relevant Boltzmann equations are Eqs. (|5.22I) - (I5.24I) or 
Eqs. (|5.19p ~ (|5.2ip . depending on whether spectator processes are taken into account 
or not. The relevance of flavour effects above T = 10 12 GeV is a novel effect specific 
to scalar triplet leptogenesis, which has been overlooked in previous works. 

• in the regime 10 9 GeV < T < 10 12 GeV, where the tau Yukawa coupling is in equi¬ 
librium while the muon and electron Yukawa couplings are not, the flavour-covariant 
computation involves the 2x2 density matrix A 9 ^ ( a , (3 = e, fi) and the asymmetry 
A r , whose evolutions are governed by the Boltzmann equations (15.51) and (15.61) . This 
computation is compared with the 2 -flavour approximation in which the 2 x 2 density 
matrix A 9 ^ is replaced by Ao (the overall asymmetry in the charges B/3 — L e and 
B/3 — L^) and the Boltzmann equations are given by Eqs. (15.271) - ([5.291) . Again, the 
relevance of flavour effects in the (£ e , £„) subspace in the temperature range where 
the electron and muon Yukawa coupligs are out of equilibrium is a new effect specific 
to scalar triplet leptogenesis. 

Below T = 10 9 GeV, flavour covariance is completely broken by fast processes induced by 
the tau and muon Yukawa couplings and the relevant dynamical parameters are A e , A^ 
and A t , whose evolution is governed by the Boltzmann equations (15.81) (15.101) . 


6.1 Parameters 

The parameters of scalar triplet leptogenesis are the triplet mass M A and its couplings 
to Higgs (n or \h = \n\ /Ma) and lepton doublets (f a p), as well as the coefficients of 
the effective dimension-5 operators (12.91) . n a p. These parameters are not all independent, 
however, as they are related to the neutrino mass matrix by 

mi + mu = rn v , ( 6 . 1 ) 

where m A and mu are defined by Eqs. (12.41) and (12.101) . respectively, and m u = U*D V U' 
in the charged lepton mass eigenstate basis, with D v = Diag(mi, m^, m 3 ), m* (i = 1,2,3) 
the neutrino masses and U the PMNS matrix. Thus, for given neutrino parameters, the 
Ka ^s are fully determined once Ma, A h and the f a p s are fixed. It is then convenient to 
write the triplet contribution to the neutrino mass matrix as 

m A = U*V*D A V ] U\ (6.2) 


where D A is a diagonal matrix with real positive entries and V is the unitary matrix that 
relates the lepton doublet basis in which m u is diagonal to the one in which m A = D A . 
We adopt the following parametrization for V: 


le 0 0 \ le iai 0 0 \ 

V = | 0 e~^ 2 0 23, ^>13, 012 , 7 ) 0 1 0 j , 

\ 0 0 e-ifej \ 0 0 e ia2 J 


U' 


( C 12 C 13 
—S12C23 — C12S13S23 
S12S23 — C12S13C23 e* 7 


S 12 C 13 

C12C23 — S12S13S23 e* 7 
—C12S23 — S12S13C23 e* 7 


513 \ 
Cl3S23e* 7 , 

Cl3C23e V 


(6.3) 


(6.4) 
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in which c ?J - = cos faj and = sin Note that the phases 01,2,3 are physical and 
cannot be removed by rephasing the lepton doublets, since this freedom has already been 
used to write the neutrino mass matrix in Eq. (16.11) in the standard phase convention. 

The flavour structure of the couplings f a p (or equivalently, of the matrix m/s) plays 
a crucial role in leptogenesis as it determines the flavour-dependent CP asymmetries and 
washout rates. Given the large number of parameters involved, we shall consider suitably 
chosen ansatze for m/s in the following. Before presenting them, let us first discuss on a 
qualitative basis how flavour effects depend on the triplet couplings. 


6.2 Qualitative discussion of flavour effects 

In the flavour-covariant formalism used in this paper, the violation of CP in triplet decays 
is described by an hermitian CP-asymmetry matrix £ a p given by Eq. (13.581) . Using the 
fact that m v = m/s + mu, one can rewrite it as 

o7TZ V z m a 


where m/s has been defined previously as m/s = y tr(m^mA) • Using the Cauchy-Schwarz 
theorem, one can derive an upper bound on the total CP asymmetry e/s = tr£ [[26]: 


£A 


< 


1 Ma 

47T V 2 


\j B/Bh rrijy , 


( 6 . 6 ) 


where fh u = y tr (mtm u ). This bound is saturated for ttt-a = iCm v with C real: the total 
CP asymmetry in triplet decays is maximal when the f a p couplings and the neutrino mass 
matrix entries have the same flavour structure and an overall phase difference ±7t/2. For 
C > 0, this gives 


£a8 — 


47 r v 


Ma / d ~ {m u ml) a p 

—y V -=- 

m v 


ca = - 


1 M a 

47T V 2 


V B/B h 


m v 


(6.7) 


In addition to a large enough CP asymmetry, successful scalar triplet leptogenesis 
requires that decays dominate over annihilations around z 1 (T r\_/ Ma), hence the 
triplet decays are typically fast. In spite of this, a large lepton asymmetry can develop 
when some decay channels are slower than the expansion of the universe. For instance, 
in the single flavour approximation, it was shown in Ref. [26] that a large efficiency is 
obtained when one of the two decay modes A —>• HH or A —> ti is out of equilibrium, 
i.e. when A u A^ or A 1 <C A h- As we are going to see, when lepton flavour effects are 
taken into account, it is also possible to reach an order one efficiency in the case Xe ~ A h, 
provided that some flavoured decay channel A —>• I a ip is slower than the expansion of 
the universe. By contrast, when A / 3> A h, all decay channels into antileptons are fast 
and only the decays into Higgs bosons can occur out of equilibrium (barring a strong 
hierarchy among the couplings f a p). In this case, asymmetric triplet and antitriplet 
decays generate an asymmetry in the Higgs sector in the first place, which is accompanied 
by a B — L asymmetry of comparable size when the triplet abundance drops [26| (see 
e.g. Eq. (14.161) . which holds above T = 10 1 ' 5 GeV, and analogous relations valid in other 


















temperature regimes in Section E|) . Lepton flavour does not play a prominent role in this 
mechanism, hence no significant difference between the single flavour approximation and 
the flavoured computation is expected in this casJ^I. A similar conclusion can be reached 
in the opposite case Xe <C A h- Indeed, for A^ small enough, the washout of the flavoured 
lepton asymmetries can be neglected and the Boltzmann equation for A Q( g simplifies to 

SHZ ~^F ~ ~ (i* “ ( 6 - 8 ) 

Taking the trace of this equation, one recovers the single flavour approximation (with 
A = tr(A a/ g) = A b-l)' 

sffz ^ =: "(J§“ 1 ) 7DeA ' (6,9) 

Hence, in the case A^ <C A h too, no significant effect of lepton flavour is expected. 


6.3 Ansatze and inputs 

The above discussion suggests that lepton flavour effects will not play a prominent role 
when the scalar triplet gives a subdominant contribution to neutrino masses (i.e. ffiA <C 
rfi u ), because in this case either T(A —>• HH) or T(A —>• ££) is smaller than the expansion 
rate of the universal so that a large efficiency is already reached in the single flavour 
approximation. For this reason, we will only consider ansatze satisfying the condition 
mi = m u . Specifically, we will study the following ansatze (explicit expressions for m a 
are given in the charged lepton mass eigenstate basis): 

• ansatz 1: ttt-a = im u 

/ y/l — x 2 y 0 

• ansatz 2: toa = irh u U* 0 xy 

V o o 

where x and y can vary between 0 and 1. Ansatz 1 is a particular case of Ansatz 2 
corresponding to xq = m 2 / y/m\ + m\ and yo = m\ + m^/rhy . 

• ansatz 3 : ttt-a = U*V* DaV^U^ with ^23 = ^>13 = 0 and «i = «2 + 7 = P3 = 0, i.e. 

(e~^ 0 0 \ / cia a 12 0 \ 

V = I 0 e~^ 2 0 -aia C 12 0 , (6.10) 

\ 0 0 1) \ 0 0 1/ 

with two different options for Da and the phases /3\, fo'- 

14 An obvious exception to this general statement is when the total CP asymmetry in triplet decays vanishes, 
while the flavour-dependent CP asymmetries do not (or more generally when ca is small due to strong 

cancellations between the diagonal entries of £). In particular, it does not apply to the scenario of purely 
flavoured leptogenesis discussed in Refs. uni ei], which includes traceless contributions to £ a p that are not 
considered in this paper (see the comment in Section ]2|). 

15 One can see this by noting that K^Kh ~ 550 (toa/0.05 eV) 2 , where K a = T(A —> aa)/H{MA) and the 
Hubble rate is given by H(T ) = 1.66 yfg££T 2 /Mp with g » = g* M = 106.75. The case in which both Ke < 1 and 
Kh < 1 is not interesting for leptogenesis, because triplets and antitriplets would annihilate before decaying. 
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ansatz 3-a: Z?a = Dv 


mi 0 0 

0 m 2 0 
0 0 m 3 


, P 2 = 0, while f3\ and <f> 12 can 


vary; 

ansatz 3-b: D\ 


xrriQ 0 0 \ 

0 y /1 — x 2 m 0 0 I , /3i = (3 2 = vr/4 , 

0 0 m3 J 


with mo = a/ m 2 + m^ , while x and (j) 12 can vary (0 < x < 1). Note that 
-Da = -Djy for x = mi /a/ m 2 + m| ■ 


Ansatz 1 depends only on 2 parameters, which can be chosen to be Ma and A h- The 
other ansatze involve 2 additional parameters that control the flavour structure of the 
couplings f al 3 : x and y for Ansatz 2, c/>i 2 and /3\ for Ansatz 3-a, cj) 12 and x for Ansatz 3-b. 

We also need to specify the neutrino mass parameters, which serve as inputs in the 
above ansatze and determine the couplings through Eq. (jG.ljh For definiteness, we 
assume a normal hierarchy. The running of the neutrino mass matrix up to the triplet 
mass scale Ma is performed using the Mathematica package REAP [53]. For a hierarchical 
neutrino spectrum, the only effect of running is to multiply the eigenvalues of the neutrino 
mass matrix by a common factor r [53], which for Ma around 10 9 -10 12 GeV lies between 
1.2 and 1.4. We therefore define the squared mass differences as AmJi = r 2 Am 2 un and 
Am§! = r 2 Arn 2 tm , where Am 2 un and Am 2 tm are the parameters extracted from oscillation 
experiments, and set the lightest eigenvalue mi = 10 ~ 3 eV at the scale Ma (which corre¬ 
sponds to a lightest neutrino mass in the range (0.7-0.8) x 10 ~ 3 eV). For the oscillations 
parameters, we take Am 2 un = 7.59 x 10 -5 eV 2 , Am 2 tm = 2.47 x 10 -3 eV 2 , sin 2 612 = 0.30, 
sin 2 623 = 0.42 and sin 2 @13 = 0.023, within ler of the most recent fit m- Finally, we set 
all CP-violating phases of the PMNS matrix to zero. 

Let us discuss these ansatze in turn. Ansatz 1 maximizes the total CP asymmetry in 
triplet decays, but due to the hierarchy among the couplings f a p implied by the relation 
m-A = im u , it is not expected to lead to very strong flavour effects. Indeed, taking into 
account the hierarchy mi, m 2 <C m 3 , one can write, in the neutrino mass eigenstate basis, 


(it 0 o\ (0 0 o\ (0 0 0 \ 

f = iXi 0 ^ 0 « *A* 0 0 0 , £ ~ 0 0 0 . (6.11) 

\o 0 m.) \o 0 l) \o 0 e A J 

Therefore, as a first approximation, leptogenesis can be described in terms of a single 
flavour (the tau flavour). But since the hierarchy between neutrino masses is not very 
strong, the difference between the single flavour approximation and the result of the 
flavoured computation can still be significant, as we are going to see. 

In Ansatz 2, the couplings / Q( g are diagonal in the neutrino mass eigenstate basis as in 
Ansatz 1 , but the diagonal entries f aa are no longer proportional to neutrino masses. As 
a result ca is not maximal, but since f a p and £ a p do not necessarily have the hierarchical 
structure of Eq. (16.lip , flavour effects can be more important than in Ansatz 1 and 
eventually lead to a larger baryon asymmetry. 

Ansatz 3 allows to investigate the effect of a misalignment between the neutrino 
mass matrix m u and the triplet couplings to leptons f a p, i.e. of a relative rotation be¬ 
tween the lepton doublet basis in which m v = D u and the one in wich the f a p s are 
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Figure 5: Comoving number density of triplets and antitriplets £a and asymmetries Aa, A h 
and Ab-l = Tr (A a( g) as a function of z — M&/T, for Ma = 5 x 10 12 GeV, A h = 0.1 and 
7tia = im v . Also shown is the departure of £a from its equilibrium value. Asymmetries are 
plotted in units of the total CP asymmetry in triplet decays ca- 

diagonal. We consider two subcases, Da = D u = Diag(mi, m 2 ,m 3 ) (Ansatz 3-a) and 
Da = Diag(xmo, Vl — x 2 mo, m 3 ), with mo = y/m‘f + m^ (Ansatz 3-b). The source 
terms are concentrated in the (A,£ 2 ) subspace (i.e. £13 = £23 = £33 = 0 in the neutrino 
mass eigenstate basis), in which the washout is weaker since \f a p\ I/ 33 I for a, ^ = 1,2 
and /13 = /23 = 0. In this case, the unflavoured computation (in which the washout is 
controlled by A £ ~ I/ 33 I) is likely to be a very bad approximation to the full, flavoured 
computation. 

6.4 Numerical results 

After solving numerically the appropriate Boltzmann equations with initial equilibrium 
abundances for the triplet^], leptons and Higgs bosons, one is left with the final value of 
the asymmetries A a = A B/ 3 -L a - The baryon asymmetry is then given by [56] 

19 10 

Ab = —A B -l = (6.12) 

a 

where we have assumed that sphalerons go out of equilibrium below the electroweak 
phase transition, as suggested by recent lattice computations m- One can also express 
the baryon asymmetry of the universe in terms of the more familiar baryon-to-photon 
ratio t]b = n^/n 7 = 7.04 Ab, whose observed value from recent Planck data is [58] 

^ bs - = (6.10 ± 0.06) x 10 " 10 ( 68 % C.L.). (6.13) 

16 This assumption is justified by the fact that gauge interactions bring the population of triplets and an- 
titriplets into equilibrium very quickly, except for very large triplet masses. 
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Figure 6: Isocurves of the baryon-to-photon ratio n^/n 7 in the ( x,y) plane for Ansatz 2. Left 
panel: Ma = 5 x 10 12 GeV, A h = 0.4. Right panel: Ma = 10 11 GeV, A h = 0.025. 


As can be seen from Fig.0 the B—L asymmetry is already stabilized around T ~ Ma/100. 

We are now ready to present our numerical results. To reduce the number of free 
parameters, we concentrate on the ansatze defined in the previous subsection and fix 
the neutrino parameters as specifed there. Let us first consider Ansatz 2, which unlike 
Ansatz 1 does not maximize ea, but allows for larger flavour effects. Since this ansatz 
gives same-sign source terms in the Boltzmann equations, flavour effects are more likely to 
be relevant when Bi ~ Bh', hence, for a given triplet mass, we choose the parameter A h 
such that Xf ~ A u- Then we study the dependence of the baryon-to-photon ratio ns/n 7 
on x and y, the parameters that control the relative sizes of the couplings / Q( g. The result 
of the full computation (i.e. the computation involving the flavour-covariant Boltzmann 
equations with spectator processes included) is shown in Fig. [6] for two different values 
of the triplet mass, Ma = 5 x 10 12 GeV and Ma = 10 11 GeV. For Ma = 5 x 10 12 GeV, 
ns/n 7 exceeds the observed value everywhere, but this can by cured by choosing a smaller 
overall phase in tua- We can see that the final baryon asymmetry is maximal for x 
small and y close to 1, more precisely around ( x,y ) = (0.05,0.95), corresponding to 
Da = Diag (0.949, 0.0475, 0.312) . The baryon-to-photon ratio in this region is typically 

enhanced by a factor 2 or 3 with respect to Ansatz 1, although the total CP asymmetry 
is smaller. These conclusions are relatively stable under variations of the triplet mass as 
long as Bi Bh- 

We now proceed to study Ansatz 3, in which the CP asymmetry is concentrated in 
the (£ 1 ,^ 2 ) subspace (i.e. £13 = £23 = £33 = 0 in the neutrino mass eigenstate basis). 
Let us first consider Ansatz 3-a, in which Da = D v while (3\ and ^>12 can vary, all other 
phases and mixing angles in the unitary matrix V being zero. We choose a triplet mass 
in the range [10 9 GeV, 10 12 GeV], and compare the result of the full computation with 
the one of the 2-flavour approximation, as described at the beginning of this section. 
As can be seen from Fig. 0 the baryon-to-photon ratio computed in a flavour-covariant 
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Figure 7: Left panel: the baryon-to-photon ratio ne/n 7 predicted by the full computation (red 
dashed lines) and by the 2-flavour approximation (blue dotted lines) as a function of j3i and 
0i2 for Ansatz 3-a. The coloured areas indicate where the observed baryon asymmetry can be 
reproduced. Right panel: ratio of the full computation over the 2-flavour calculation of the 
baryon asymmetry. Along the solid curve, n b /n 7 vanishes in the 2-flavour approximation. In 
both panels, Ma = 5 x 10 11 GeV and A h = 0.1. 


way can be enhanced by up to a factor of several hundreds with respect to the 2-flavour 
approximation. Indeed, important flavour effects in the e—fi sector are missed in this 
approximation. This can be understood by examining a specific point of the parameter 
space, for instance 0i2 = vr/6 and /3\ = tt/2. In the full computation, the sources terms 
in the Boltzmann equations (15.51) and (|5.6I) are proportional to 


7-10.36 -2.496A 
V—2.496 6.846 ) X 0 


Err = 3.513 x 10“ 8 , 


(6.14) 


where for definiteness we have written £° in the (£ e , 1^) basis, while the flavour dependence 
of the washout is controlled by the triplet couplings to leptons: 


/l/ee| 

1 /e/,| 

1 fer | \ 

/ 1.428 

3.477 

1.373A 


( \fpe\ 


\Ur\ - 

= 3.477 

9.288 

9.394 x 10 -3 . 

(6.15) 

Vl/rel 

I/t/,1 

\frr\J 

\1.373 

9.394 

13.36/ 



Thus we have a large source term for a flavour asymmetry (namely in the e flavour) that 
is only weakly washed out by inverse decays. As a result, a large baryon asymmetry is 
generated ( ub/u 7 = 7.277x 10~ 10 ) even though the total CP asymmetry vanishes (indeed, 
for the particular point considered ca = tr£° + S TT = 0). By contrast, in the 2-flavour 
approximation, the CP asymmetries appearing in the Boltzmann equations (15.281) are 

e 0 = -3.513 x 1(T 8 = -e T , (6.16) 
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Figure 8 : Left panel: the baryon-to-photon ratio Ub/ti 7 predicted by the full computation (red 
dashed lines) and by the 2 -flavour approximation (blue dotted lines) as a function of x and 
0 i 2 for Ansatz 3-b. The coloured areas indicate where the observed baryon asymmetry can be 
reproduced. Right panel: ratio of the full computation over the 2-flavour calculation of the 
baryon asymmetry. In both panels, Ma = 5 x 10 11 GeV and A h — 0.1. 


and the washout is controlled by 

(foo for\ = ( 10.69 9.567A 

V/0t \frr\J V 9 ' 567 13 - 36 / 


(6.17) 


Therefore, the asymmetries stored in the two flavours are washed out with a comparable 
strength, and the resulting baryon asymmetry (ns/n 7 = 1.617 x 10~ 12 ) is much smaller 
than in the full flavoured computation. 

A similar (but numerically not as strong) enhancement of the full computation result 
with respect to the 2-flavour approximation can be observed for Ansatz 3-b, in which 
Da = Diag(x mo, VT — x 2 mo, m 3 ) with mo = y m\ + m\ and the phases j3\ and /?2 are 
set to 7 t/ 4. This is shown in Fig. [HI 

In Figs. l9l and 1101 we compare the relative impacts of spectator processes and flavour 
covariance on the generated baryon asymmetry in the temperature regimes T > 10 12 GeV 
and 10 9 GeV < T < 10 12 GeV, respectively. We consider two ansatze designed to produce 
a large baryon asymmetry (even exceeding the observed value when the CP-violating 
phases are chosen to be large, as is the case here), namely Ansatz 1, which maximizes the 
total CP asymmetry, and Ansatz 2 with (x,y) = (0.05,0.95), a parameter choice that has 
been shown to maximize the final baryon asymmetry in Fig. [ 6 j In Fig. [9l the triplet mass 
has been chosen to be Ma = 5x 10 12 GeV, so that most of the B—L asymmetry is produced 
at T > 10 12 GeV. For Ansatz 1 (i.e. toa = ifn„), flavour effects and spectator processes 
have a quantitatively similar impact on the baryon-to-photon ratio, while for Ansatz 2 
flavour effects strongly dominate over spectator processes such as QCD sphalerons or top 
quark Yukawa interactions (especially for 1CT 2 < < 0.5). In this case, neglecting 
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Figure 9: Baryon-to-photon ratio ns/n 7 as a function of A e for Ma = 5 x 10 12 GeV, assuming 
Ansatz 1 (left panel) or Ansatz 2 with (x,y) = (0.05,0.95) (right panel). The red lines indicate 
the result of the flavour-covariant computation involving the 3x3 density matrix A Q( g, with 
(solid red line) or without (dashed-dotted red line) spectator processes taken into account, 
whereas the blue lines indicate the result of the single flavour approximation, taking spectator 
processes into account (blue dashed line) or not (blue dotted line). The equality of branching 
ratios Be = Bh is realized for A e ~ 0.15. 




Figure 10: Baryon-to-photon ratio ub/ti 7 as a function of A e for Ma = 10 11 GeV, assuming 
Ansatz 1 (left panel) or Ansatz 2 with (x,y) = (0.05,0.95) (right panel). The red lines indicate 
the result of the flavour-covariant computation involving the 2x2 density matrix A^, with 
(solid red line) or without (dashed-dotted red line) spectator processes taken into account, 
whereas the blue lines indicate the result of the 2-flavour approximation, taking spectator 
processes into account (blue dashed line) or not (blue dotted line). The equality of branching 
ratios Be = Bh is realized for Xe ~ 0.021. 
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Figure 11: Isocurves of the baryon-to-photon ratio %/n-y in the (A^,Ma) plane obtained 
performing the full computation, assuming Ansatz 1 (left panel) or Ansatz 2 with (x, y) = 
(0.05, 0.95) (right panel). The coloured regions indicate where the observed baryon asymmetry 
can be reproduced in the full computation (light red shading) or in the single flavour approxi¬ 
mation with spectator processes neglected (dark blue shading). The solid black line corresponds 
to Bi = Bh■ Also shown are the regions where A h is greater than 1 or 47 t. 

spectator processes and including flavour effects gives a much more accurate result than 
doing the opposite. In Fig. [TOl the triplet mass is Ma = 10 11 GeV, hence the B — L 
asymmetry is essentially generated in the temperature regime 10 9 GeV < T < 10 12 GeV 
where the tau Yukawa coupling is in equilibrium, but the muon Yukawa coupling is not. 

For Ansatz 1, the 2-flavour calculation turns out to be a rather good approximation to 
the flavour-covariant computation, while neglecting spectator processes gives a very bad 
estimate, except for small or large values of A £. For Ansatz 2, on the contrary, both flavour 
covariance and spectator processes have a significant impact on the baryon-to-photon ratio 
(except again for extreme values of A^), and neglecting one of them underestimates the 
result by up to a factor 2. The 2-flavour approximation without spectator processes 
actually gives a much larger disagreement with the full computation. 

Finally, we study in Fig. [lT]the dependence of the generated baryon asymmetry on A^ 
and Ma, both for Ansatz 1 and for Ansatz 2 with ( x,y ) = (0.05,0.95). The computation 
is performed assuming that the third generation Yukawa couplings as well as the charm 
Yukawa coupling are in equilibrium, which strictly speaking is true only in the temperature 
range 10 9 GeV < T < 10 12 GeV. From Fig. [TT] one can conclude that successful scalar 
triplet leptogenesis is possible for a triplet mass as low as 4.4 x 10 10 GeV, to be compared 
with 1.2 x 10 11 GeV in the approximation where flavour effects and spectator processes 
are neglected. Other assumptions about the flavour structure of the triplet couplings 
to leptons may allow for a lighter scalar triplet. For comparison, we quote the lower 
bounds found by Ref. [26] in the single flavour approximation with spectator processes 
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neglected: Ma > 1.3 x 10 11 GeV in the case rliA = 0.05 eV (« m„ for a hierarchical 
neutrino mass spectrum), in agreement with our result, and Ma > 2.8 x 10 10 GeV in the 
case ffiA = 0.001 eV <C fh v . Although we did not consider ansatze satisfying toa rh u , 
because flavour effects are comparatively smaller in this case, we also expect the latter 
bound to be lowered by the inclusion of flavour and spectator processes. 

7 Conclusions 

In this paper, we have shown how to consistently include the effects of the different 
lepton flavours in scalar triplet leptogenesis. When charged lepton Yukawa interactions 
are out of equilibrium at the time of leptogenesis, i.e. when the lepton asymmetry is 
generated at high temperature, the proper treatment of flavour effects involves a 3 x 3 
density matrix, whose evolution is governed by flavour-covariant Boltzmann equations. 
The often-used single flavour approximation, which gives rather accurate results in the 
standard leptogenesis scenario with right-handed neutrinos, leads to predictions for the 
generated baryon asymmetry that can depart by a large amount from the flavour-covariant 
computation. In the intermediate temperature regime where the tau Yukawa coupling is 
in equilibrium but the muon and the electron ones are not, the 3x3 density matrix can 
be replaced by the asymmetry stored in the tau lepton doublet and by a 2 x 2 density 
matrix describing the flavour asymmetries in the (£ e , £n) subspace and their quantum 
correlations. In this case too, the 2-flavour calculation in which the (7 e , 7 /t ) subspace is 
described by a single flavour £q does not in general give a good approximation of the 
flavour-covariant computation. Finally, when the tau and muon Yukawa couplings are in 
equilibrium, flavour covariance is completely broken and the dynamics of leptogenesis is 
described in terms of the asymmetries stored in the three lepton doublets £ e , £ il and i T . 

We performed a numerical study of the impact of flavour effects and spectator pro¬ 
cesses on the generated baryon asymmetry for judiciously chosen ansatze, and compared 
the flavour-covariant computation with flavour non-covariant approximations used in the 
literature, with or without spectator processes included. We found discrepancies in the 
predictions for the baryon asymmetry ranging from an order one factor to two orders of 
magnitude. In particular, we showed that successful leptogenesis can easily be achieved 
when the decays of the triplet into leptons and Higgs bosons occur at a similar rate, while 
it would require a significantly heavier triplet in the single flavour approximation. As 
a result, the minimal triplet mass allowed by successful leptogenesis is lowered by the 
inclusion of flavour effects, from 1.2 x 10 11 GeV to 4.4 x 10 10 GeV in the case where the 
scalar triplet and the additional heavy states give comparable contributions to neutrino 
masses. 

Throughout this paper, we worked in a framework in which the contribution of the 
additional heavy states to the CP asymmetries is parametrized by an effective dimension-5 
operator, but the procedure we used to derive the flavour-covariant Boltzmann equations 
can also be applied to explicit models with e.g. several scalar triplets, or with a scalar 
triplet and right-handed neutrinos. The formalism employed in this paper can also be 
used to study the scenario of purely flavoured leptogenesis [Sol |3l] in a flavour-covariant 
way. In the framework used in this paper, this would require the addition of the effective 
four-lepton operator (|2.11h . 
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A Reaction densities 


We summarize in this appendix some useful formulae for the reaction densities used in this 
paper. Let us first recall the general expression for the (thermally averaged) space-time 
density of a general reaction a + b + ...—ti+j + .... Neglecting Bose enhancement and 
Pauli blocking factors: 


7 (a + b + ... ->• i + j + ...) 


d 3 Pa eq _d 3 pb_ eq d 3 p.j d 3 pj 

(27r)32 Wp - a Pa 2u fb (2tt)3 Pb ''' 2u Pi (2vr)3 2w ft . (2 tt)3 

|A4| 2 (2ir) 4 5^\p a +Pb + ■ ■ ■ - Pi ~ Pj ~ ■ • •), (A.l) 


where |A4| 2 is the squared matrix element summed over the internal degrees of freedom 
of the initial and final states, and (p) is the phase-space distribution function of the 
particle X at kinetic and chemical equilibrium, which only depends on the bosonic or 
fermionic nature of X: 

= (A ' 2) 

In the following all computations are done using the Maxwell-Boltzmann statistics, i.e. 
neglecting the difference between bosons and fermions: 

4’W = Pb(p) = ‘•-' 3(£_ " ) = Pm’bW • (A.3) 


With this approximation, the space-time density of triplet and antitriplet decays can be 
written as 

1D = K 2 (z) Ta ’ ^ A ' 4 ^ 

where Ea = (n^ + n^)/s is the comoving number density of triplets and antitriplets, 
K\ 2 (z) are modified Bessel functions of the second kind, z = M^/T and Ta is the triplet 
decay width. For 2 —>■ 2 scatterings, one has 

7 (a + b^i + j) = / dss 1/2 a(s)K 1 (^Y) , (A.5) 

where a(s) = 2s \(l,m 2 /s,ml/s) cr(s) is the reduced cross-section summed over the in¬ 
ternal degrees of freedom of initial and final particles. 

When computing the space-time density of a 2 —>• 2 scattering, one must take care 
to properly subtract the contribution of on-shell intermediate particles, which is already 
taken into account in decays and inverse decays [59]. When the resonance occurs in the 
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s-channel, one can compute the subtracted reaction density by taking away the resonant 
part from the squared propagator in the narrow-width approximation EM: 


\D\ 


\D \ 2 - 


7T 

Mr 


5(s- M 2 ) , 


(A. 6 ) 


where M and T are the mass and width of the intermediate particle, and D is the propa¬ 
gator of the intermediate state in the Breit-Wigner form: 

1 


D = -- . (A.7) 

s — M 2 + iMT v ' 

We have to do this for the computation of ^{tt —>• HH) and 7 (££ —> ££), which receive 
a contribution from s-channel triplet exchange. Similarly, when computing y(£A —>• ^A), 
one has to subtract the contribution of real intermediate leptons in the u-charmel, corre¬ 
sponding to the two processes A —>• ££ and it —» A. Following Ref. [5], we perform the 
subtraction in the following way: 


\D\ 


\D\ 2 - 


TT 


ET 


th 


-5 (u) 


(A. 8 ) 


where E is the energy of the intermediate lepton, r t h its thermal width, and the propagator 
of the intermediate lepton is 

D = - . (A.9) 


u + iE r th 

As in Ref. [5], we use the following representation of 5{x) in numerical computations: 


5{x) = 


(A.10) 


7 r (x 2 + e 2 ) 2 ’ 

where e is a small number. In the limit of small width, which we assume to be valid, one 
can simply set e = T/M for a resonance in the s-channel with an intermediate particle of 
mass M and width T, and £T t h = eM| for a resonance in the u-channel, with any small 
value for e. With this choice, one can compute the subtracted space-time densities for the 
various 2 —2 scatterings. 

In what follows, we note x = s/M The reduced cross-sections for the scatterings 
involving 2 leptons and 2 Higgs bosons are, after subtracting the contribution of the 
on-shell intermediate triplet, 

(1 — x) 2 — e 2 
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for ti —>• HH, and 
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(A.12) 
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for IH —>• IH, where the superscripts A, 'H and X refer to the contributions of the scalar 
triplet, of the D = 5 operator (12.91) responsible for mu and to the interference between 
the two contributions, respectively. For the 2 lepton-2 lepton scatterings ££ —>• li and 
il —>• £&, we obtain 


and 


cr, 


M 


3 \ 4 2 (1 — x) 2 — e 2 
32vr 1 [(! - x f + e 2] 2 


(A.13) 






2 

1 + x 


4 \ 

— ln(l + x) . 
x J 


(A.14) 


We did not include the gauge contributions since they do not violate flavour, hence they 
drop from the Boltzmann equations. Finally, the reduced cross-sections for the scatterings 
(A —>• £A, £A —>• iA and £i —> A A are given by, after subtracting the contribution of the 
on-shell intermediate lepton in the ti-channel, 
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yj{x — 4)x \x — 2 — y/x(x — 4) 


(A.15) 
(A.16) 

(A.17) 


With the above expressions, one can compute numerically the reaction densities 7 
74 £ and 7 ^a that appear in the washout terms and W^, respectively (see Sub¬ 

section [2T3] for the definition of these reaction densities and Eqs. (|3.52l) . (I3.53P and (I3.54|) 
for the washout terms). We also need the space-time density of triplet-antitriplet annihi¬ 
lations into Standard Model particles, 7 , 4 , which enters the Boltzmann equation for Sa, 
Eq. (13.601) . The contribution of gauge scatterings to 7,4 has been computed in Ref. [ 26] . 
There are also subleading contributions proportional to \j, A#, A 2 g 2 and A 2 H g 2 , which we 
do not include in our computation. 

Fig.ll2lshows the typical magnitude of the various reaction densities and their evolution 
as a function of z = M&/T (here for Ma = 5 x 10 12 GeV, toa = im v and A h = 0 . 2 ). 
A ratio 7 smaller than 1 indicates that the reaction R is slow on a cosmological 
time scale. In this example, most scatterings are negligible, except for the ones involving 
Higgs bosons; however, the latter become slower than the expansion of the universe at the 
time where decays start to dominate over annihilations, a necesary condition for a large 
lepton asymmetry to develop [ 26l 132] , For the present choice of parameters, however, the 
third Sakharov condition would not be satisfied in the single flavour case, as both decays 
into Higgs bosons and into leptons are in equilibrium at the time of leptogenesis (i.e. at 
z < 10, before the triplet abundance becomes strongly Boltzmann-suppressed). Taking 
into account flavour effects leads to an enhanced efficiency because some “flavoured” decay 
channels A —>■ £ a £p are significantly slower than others. 
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Figure 12 : Space-time densities of decays (Be 7 ^, Bh^d) and scatterings ( 7 , 4 , 7 ^, 7 ^ and 7 ^) 
in units of ifn 7 , as a function of z — Ma/T , for the values of the parameters indicated in the 
text. 
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